EXPLICIT METHODS FOR HILBERT MODULAR FORMS 



LASSINA DEMBELE AND JOHN VOIGHT 

Abstract. We exhibit algorithms to compute systems of Hecke eigenvalues for spaces of 
Hilbert modular forms over a totally real field. We provide many explicit examples as well 
as applications to modularity and Galois representations. 
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The study of modular forms remains a dominant theme in modern number theory, a conse- 
quence of their intrinsic appeal as well as their applications to a wide variety of mathematical 
problems. This subject has seen dramatic progress during the past half-century in an en- 
vironment where both abstract theory and explicit computation have developed in parallel. 
Experiments will remain an essential tool in the years ahead, especially as we turn from 
classical contexts to less familiar terrain. 

In this article, we discuss methods for explicitly computing spaces of Hilbert modular 
forms, refashioning algorithms over Q to the setting of totally real fields. Saving definitions 
for the sections that follow, we state our main result. 

Theorem. There exists an algorithm that, given a totally real field F, a nonzero ideal 9t of 
the ring of integers of F, and a weight k e (Z>2)t F: ^, computes the space S^OT) of Hilbert 
cusp forms of weight k and level OT over F as a Hecke module. 

This theorem is the work of the first author [15] together with Donnelly [IE] combined 
with work of the second author [69] together with Greenberg [30] . 

The outline of this article is as follows. After briefly recalling methods for classical (el- 
liptic) modular forms in §1, we introduce our results for computing Hilbert modular forms 
in the simplest case (of parallel weight 2 over a totally real field of strict class number 1) in 
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§2. In brief, our methods employ the Jacquet-Langlands correspondence to relate spaces of 
Hilbert modular forms to spaces of quaternionic modular forms that are more amenable to 
computation; we discuss this matter in §3, and consider two approaches (definite and indefi- 
nite) in §§4-5. In §6 we consider several examples in detail. Having thereby established the 
main ideas, we turn to an adelic description of Hilbert modular forms and their quaternionic 
incarnations in §7, then give a complete and general description of our algorithms in adelic 
language in §§8-9. 

Although it is our intention to keep these notes as self-contained as possible, we will as- 
sume that the reader has a basic familiarity with classical modular forms and the methods 
employed to compute them. The algorithms exhibited below have been implemented in the 
computer algebra system Magma [3] and our examples are computed using this implementa- 
tion. Donnelly and the second author [21] are using this implementation to compute Hecke 
data for thousands of forms over totally real fields up to degree 6. 

These notes arose from lectures at the Centre de Recerca Matematica (CRM) in Barcelona; 
it is our great pleasure to thank the CRM for the invitation to speak and the hospitality 
of the organizers, Luis Dieulefait and Victor Rotger. The authors would also like to thank 
Matthew Greenberg, Ariel Pacetti, Aurel Page, Jeroen Sijsling, and the referee for many 
helpful comments as well as Benedict Gross for his remarks which we include at the end of 
Section 6. The first author is supported by a Marie-Curie Fellowship, and the second author 
by an NSF Grant No. DMS-0901971. 



I. Classical (elliptic) modular forms 

To motivate our study of Hilbert modular forms, we begin by briefly considering algorithms 
for classical (elliptic) modular forms. For a more detailed introduction to modular forms, 
see the books by Darmon [llj and Diamond and Shurman [19], and for more information 
on computational aspects see Cremona (9|, Kilford |^D] . Stein j6l], and the many references 
contained therein. 

Let U = {x + yi G C : y > 0} denote the upper-half plane and let W = H U P X (Q) denote 
the completed upper half-plane with the cusps P X (Q). The group 

GL+(Q) = { 7 GGL 2 (Q):det 7 >0} 

acts on W by linear fractional transformations. For N G Z>o, we consider the subgroup of 
those integral matrices of determinant 1 that are upper-triangular modulo N, 

r {N) = 1 7 = r b Ae SL 2 (Z) : N | cj C GL+(Z) = SL 2 (Z) C GL+(Q). 

The group Pr (iV) = r (A^)/{±l} is a discrete subgroup of PSL 2 (IR). A modular form of 
weight k G Z >0 and level N is a holomorphic function / : Ti — > C such that 

< L1 > /(7*) = /(^) + 

for all 7 G T (N) and such that f(z) tends to a finite limit as z tends to any cusp (i.e., / is 
holomorphic at the cusps). 

One can equivalently write this as follows. For 7 G GL 2 (IR) and z G H we define j(7, z) = 
cz + d. We then define a weight k action of GL^ (Q) on the space of complex- valued functions 
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on % by 

(1.2) (f\kl)(z)= ide ^ )k ' 1 f(iz). 



Then (1.1) is equivalent to /|fc7 = / for all 7 G T (N). 

Note that the determinant factor (det7) fc_1 in our definition is different from the usual 
(det7) fe//2 , which is an analytic normalization. Consequently, the central subgroup Q x C 
GL+(Q) acts by f\ kl = 1 k ^f. 

Remark 1.3. For simplicity we treat only the case of r (A^)-level structure in this article. 
If desired, one could without difficulty extend our methods to more general level structures 
with characters, and so on. 

The C-vector space of modular forms of weight k and level N is finite-dimensional and is 
denoted M k (N). 

If / G Mjfc(iV), then f(z + 1) = f(z) so / has a Fourier expansion 

00 

(1.4) f(z) = ^ a nQ n = a o + a x q + a 2 q 2 + a 3 q 3 + ... 

n=0 

where a n G C and q = exp(2niz). 

We say that / is a cusp form if f(z) — > as z tends to any cusp (i.e., / vanishes at the 
cusps). The C-vector space of cusp forms of weight k and level N is denoted S k {N). We 
have Mk(N) = Sk(N) © E k (N) where E k (N) is spanned by the Eisenstein series of level N. 



Note that when k > 2, then (1.1) is equivalent to 

f(jz) {d{ 1Z )) k - 1 = f(z) (dzf- 1 

so one may equivalently think of a cusp form / G Sk{N) as a holomorphic differential 
(k — l)-form on the modular curve X (N) = T (N)\H* . (Because of our normalization, such 
differential forms will be global sections of the line bundle corresponding to the algebraic 
local system Sym fc ~ 2 (C 2 ), corresponding to bivariate homogeneous polynomials of degree 
k — 2. Some authors use an analytic normalization instead.) 

The spaces Mk(N) and Sj~(N) are equipped with an action of pairwise commuting diago- 
nalizable Hecke operators T n for each integer n G Z>o- The Hecke operators can be thought 
of in several different ways: they arise from correspondences on the modular curve X (N), 
as "averaging" operators over lattices of index n, or more formally from double coset de- 
compositions for the group T (N) inside SL 2 (Z). The action of the Hecke operator T n is 
determined by the action of T p for p | n, and the latter for p \ N in weight k are given simply 
by the formula 

(T p f)(z)=p^f(pz)+ 1 -Y^f( Z -^ 

(For primes p \ N one omits the first term, and there are also operators called Atkin-Lehner 
involutions.) We say therefore that Sk{N) is a Hecke module, namely, an abelian group 
equipped with an action of the Hecke algebra T = Z[T P ] P = Z[T 2 , T 3 , . . . ], a polynomial ring 
in countably many variables over Z indexed by the primes. Our Hecke modules will always 
be finite-dimensional C-vector spaces. 
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A form / G Sk{N) is an oldform (at d) if f(z) = g(dz) for some g G Sk{M) with M \ N a 
proper divisor and d \ N/M; we say / is a newform if / is a normalized eigenform which is 
orthogonal to the space of oldforms (with respect to the Petersson inner product). 

The space Sk{N) consequently has a basis of eigenforms, i.e., functions that are eigen- 
functions for each Hecke operator T n . If / is an eigenform, normalized so that a% = 1 in its 



g-expansion (1.4), then T n f = a n f. Moreover, the field Q({a„}) = E C C is a number field 
and each Hecke eigenvalue a n is an algebraic integer in E. 

In this way, the system of Hecke eigenvalues (a p ) p for a normalized eigenform / G Sk(N) 
determine the form /:%—>■ C These eigenvalues also determine the L-series 

n=l 71 V P P / P |7V V P 7 

associated to / (defined for Res > 1), as well as the £-adic Galois representations 

PU : Gal(Q/Q) GL 2 (Z £ ) 
associated to / with the property that for any prime p \ EN, we have 

Tr(p /j£ (Frob p )) = a p (f) and det(p /j£ (Frob p )) = p k ~ l . 

Several methods have been proposed for making the Hecke module Sk{N) into an object 
amenable to explicit computation. With a view to their generalizations to Hilbert modular 
forms, we mention two approaches which have seen wide application. (We neglect the method 
of graphs [17] and a method which uses the Eichler-Selberg trace formula [35].) For simplicity, 
we restrict our discussion to the case of weight k = 2. 

The first method goes by the name modular symbols and has been developed by Birch, 
Swinnerton-Dyer, Manin, Mazur, Merel, Cremona [9], Stein [M], and many others. The 
Hecke operators T p act naturally on the integral homology Hi(X (N), Z; cusps) — linear com- 
binations of paths in the completed upper half plane H* whose endpoints are cusps and whose 
images in X (N) are linear combinations of loops — and integration defines a nondegenerate 
Hecke-equivariant pairing which gives rise to an isomorphism (the Eichler-Shimura theorem) 



Hi(X (N), C; cusps) S 2 (N) © S 2 (N) 

where — denotes complex conjugation. The formalism of modular symbols then presents 
the space Hi(X (N),Z; cusps) explicitly in terms of paths in T-L* whose endpoints are cusps 
(elements of P X (Q)) and whose images in Xo(iV) are a linear combination of loops. We have 
an explicit description of the action of the Hecke operators on the space of modular symbols, 
and the Manin trick (the Euclidean algorithm) yields an algorithm for writing an arbitrary 
modular symbol as a Z-linear combination of a finite set of generating symbols, thereby 
recovering S 2 {N) Hecke module. 

The second method goes by the name Brandt matrices and goes back to Brandt, Eichler 
[231 121] j Pizer [IB], Kohel [13], and others. In this approach, a basis for S 2 (N) is obtained 
by linear combinations of theta series associated to (right) ideals in a quaternion order of 
discriminant N. These theta series are generating series which encode the number of elements 
in the ideal with a given reduced norm, and the Brandt matrices which represent the action 
of the Hecke operators are obtained via this combinatorial (counting) data. 
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2. Classical Hilbert modular forms 



We now consider the situation where the classical modular forms from the previous section 
are replaced by forms over a totally real field. References for Hilbert modular forms include 
Freitag [25J, van der Geer [26J and Goren [29|. 

Let F be a totally real field with [F : Q] = n and let Zj? be its ring of integers. The case 
n = 1 gives F = Q and this was treated in the previous section, so we assume throughout 
this section that n > 1. Let vi, . . . , v n : F — > R be the real places of F, and write Vi(x) = Xj. 
For 7 G M 2 (F) we write 7* = ^(7) G M 2 (R). 

For simplicity, in these first few sections (§§2-6) we assume that F has strict class number 
1; the general case, which is more technical, is treated in §7 and thereafter. 

The group 

GL+(F) = {7 G GL 2 (F) : det7< > for i = 1, . . . ,71} 
acts naturally on T-L n by coordinatewise linear fractional transformations 

diZi + L 



z i-)- 72; = (jiZ 
For a nonzero ideal 9T C "Lp, let 



i=l, 



r (or) = | 7 = (° ^ g gl+(z f ) : c g or} c gl+(z f ) c gl 2 (f). 

Let Pr (9T) = r (9c)/Z* C PGL+(Z F ). Then the image of YT (% under the embeddings 
7 ^ (Ti)i is a discrete subgroup of PGL^ (IR) n . 

Under the assumption that F has strict class number 1, we have 

Z* + = {x G Z* : Xi > for all z} = Z* 2 



pr„(<n) = { 7 = ( I 5 ) G SL 2 (Z F ) : c G }> /{±1} 



and hence GL^ (Z^) = Z£ SL 2 (Z^), and so alternatively we may identify 

'a 

v ,c d , 
in analogy with the case F — Q. 

Definition 2.1. ^4 Hilbert modular form of parallel weight 2 and Zei>e/ 91 is a holomorphic 
function f : H n — >■ C snc/i i/iai 

/ n tt \ £ ( a \ z \ + ^1 O-nZn + b n \ I T-r (CjZj + di) 2 \ 

2.2 / 72; = / — , . . . , — =11 —j- /(*) 

Vci^i + di c n z n + d n J det7- ' 

/or a// 7 G r (0T). 



i=i 



We denote by M 2 (9I) the space of Hilbert modular forms of parallel weight 2 and level 9t; 
it is a finite-dimensional C- vector space. The reader is warned not to confuse M 2 (9T) with 
the ring M 2 (i?) of 2 x 2-matrices over a ring R. 



Remark 2.3. There is no holomorphy condition at the cusps in Definition 2.1 as there was 
for classical modular forms. Indeed, under our assumption that [F : Q] = n > 1, this follows 
automatically from Koecher's principle [261 §1]. 
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Note also that if u G Z£, then 7(1*) = (g G GL 2 (Z,p) acts trivially on H and at the 

same time gives a vacuous condition in (2.2), explaining the appearance of the determinant 
term which was missing in the classical case. 



Analogous to (1.2), we define 

<"> (/i7)w=(n^)/w 



for / : "H™ -> C and 7 G GL+(F); then Ju} is equivalent to (/|7)(*) = /(z) for all 
7 G I'o(^). 

The group GL^"(F) also acts naturally on the cusps F 1 (F) <-» P 1 (IR) n . We say that 
/ G M 2 (01) is a cusp /orm if f(z) — > whenever z tends to a cusp, and we denote the 
space of cusp forms (of parallel weight 2 and level 01) by S^Ol). We have an orthogonal 
decomposition M 2 (01) = 5*2 (01) ffi-E 2 (01) where -E 2 (01) is spanned by Eisenstein series of level 
01; for level 01 = (1), we have dim £2(1) = #C1 + Zi?, where C1 + Z,f denotes the strict class 
group of "Lp- 

Hilbert modular forms admit Fourier expansions as follows. For a fractional ideal b of F, 

let 

b + — {x G b : Xi > for i — 1, . . . , n}. 

Let be the different of F, and let D -1 denote the inverse different. A Hilbert modular form 
/ G M 2 (01) admits a Fourier expansion 



(2.5) /(z) =a + ^ a M e 



27riTr(^2) 
M6(3- 1 )+ 

with a = if / is a cusp form. 

Let / G M 2 (01) and let n C Z^ be a nonzero ideal. Then under our hypothesis that F has 
strict class number 1, we may write n = ud^ 1 for some v G D+; we then define a n = a u . The 
transformation rule (2.2) implies that a„ does not depend on the choice of v, and we call a n 
the Fourier coefficient of / at n. 

The spaces M 2 (01) and ^2 (01) are also equipped with an action of pairwise commuting 
diagonalizable Hecke operators T n indexed by the nonzero ideals n of Zj?. For example, given 
a prime p f 01 and a totally positive generator p of p we have 

(2-6) (T p f)(z) = N(p)f(pz) + ^^/ ( Z -^) , 

where F p = Z^/p is the residue field of p; this definition is indeed independent of the choice 
of generator p. Using the notation (2.4), we can equivalently write 

(2-7) (T P f)(z)= J2 (/!*•)(*) 

where = ^ ^ and n a = ^ for a G F p . 

If / G ^(Ol) is an eigenform, normalized so that a(i) = 1, then T n f = a n f, and each 
eigenvalue a n is an algebraic integer which lies in the number field E = Q({a n }) C C (see 
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Shimura [58l Section 2]) generated by the Fourier coefficients of /. We again have notions of 
oldforms and newforms, analogously defined (so that a newform is in particular a normalized 
eigenform) . 

Associated to an eigenform / G S 2 (91) we have an L-function 

n 

and t-adic Galois representations 

p u : Gal(F/F) -> GL 2 (Z F>( ) 

for primes \ of Z F such that, for any prime p f 191, we have 

Tr^Frobp)) = a p (f) and det(p/,[(Frob p )) = iVp. 

Each of these is determined by the Hecke eigenvalues a n of /, so we are again content to 
compute S^Ol) as a Hecke module. 

We are now ready to state the first version of our main result. 

Theorem 2.8 (Dembele [15] . Greenberg-Voight [30]). There exists an algorithm which, given 
a totally real field F of strict class number 1 and a nonzero ideal 91 C Zp, computes the 
space 5*2 (OX) of Hilbert cusp forms of parallel weight 2 and level 91 over F as a Hecke module. 

In other words, there exists an explicit finite procedure which takes as input the field F 
and the ideal 91 C Z F encoded in bits (in the usual way, see e.g. Cohen [7]), and outputs a 
finite set of sequences (a p (/)) p encoding the Hecke eigenvalues for each cusp form constituent 
/ in S2 (91), where a p (/) G Ef C Q. (This algorithm will produce any finite subsequence 
in a finite amount of time, but in theory will produce the entire sequence if it is left to run 
forever.) Alternatively, this algorithm can simply output matrices for the Hecke operators 
T p ; one recovers the constituent forms using linear algebra. 

Example 2.9. Let F = Q(y/E). Then Zp = Z[w] where w = (l+y/b)/2 satisfying w 2 — w — 1 = 
0. Let 91 = (3w-14) C Z F ; we have iV(91) = 229 is prime. We compute that dim (91) = 4. 
There are 2 Hecke irreducible subspaces of dimensions 1 and 3, corresponding to newforms 
/ and g (and its Galois conjugates). We have the following table of eigenvalues; we write 
p = (p) for p G Z F . 



p 


(2) 


(w + 2) 


(3) 


(w + 3) (w- 4) 


Np 


4 


5 


9 


11 11 




-3 


-4 


-1 


-2 




t 


t 2 - At + 1 


-t 2 + 2t + 2 


t 2 - 2t - 3 -3t 2 + 8t + 1 



Here, the element t e Q satisfies t 3 — 3t 2 — 1+1 = and E = Q(t) is an S^-field of discriminant 
148. 

Recall that in the method of modular symbols, a cusp form / e S^iV) corresponds to a 
holomorphic differential (l-)form (27rz)/(z) cfe on X (iV) and so, by the theorem of Eichler- 
Shimura, arises naturally in the space H l (X (N), C). In a similar way, a Hilbert cusp form 
/ G ^(Ol) gives rise to a holomorphic differential n-form {2'Ki) n f(zi, . . . , z n ) dz\ ■ ■ ■ dz n on 
the Hilbert modular variety A (91), the desingularization of the compact space r (91)\('H n )* 
where (7i n )* = W 1 U F 1 (F). But now Xo(91) is an algebraic variety of complex dimension n 

7 



and / arises in the cohomology group if n (Xo(9T), C). Computing with higher dimensional 
varieties (and higher degree cohomology groups) is not an easy task! So we seek an alternative 
approach. 

Langlands functoriality predicts that 5 2 as a Hecke module occurs in the cohomology 
of other "modular" varieties as well. This functoriality was already evident by the fact 
that both modular symbols and their quaternionic variant, Brandt matrices, can be used 
to compute the classical space S^iV). In our situation, this functoriality is known as the 
Jacquet-Langlands correspondence, which ultimately will allow us to work with varieties of 
complex dimension 1 or by considering twisted forms of GL2 over F arising from quaternion 
algebras. In dimension 1, we will arrive at an algorithm which works in the cohomology of 
a Shimura curve, analogous to a modular curve, and thereby give a kind of analogue of 
modular symbols; in dimension 0, we generalize Brandt matrices by working with theta 
series on (totally definite) quaternion orders. 

3. Quaternionic modular forms 

In this section, we define modular forms on quaternion algebras; our main reference is 
Hida [3l]. We retain the notation of the previous section; in particular, F is a totally real 
field of degree [F : Q] = n with ring of integers TLp. 

A quaternion algebra B over F is a central simple algebra of dimension 4. Equivalently, a 
quaternion algebra B over F is an F-algebra generated by elements i, j satisfying 

(3.1) i 2 = a, j 2 = b, and ji = —ij 

for some a,b G F x ; we denote such an algebra B = For more information about 

quaternion algebras, see Vigneras [66] . 

Let B be a quaternion algebra over F. Then B has a unique involution ~~ : B — >■ B 
called conjugation such that xx G F for all x G B; we define the reduced norm of x to be 

a, 6" 



nrd(x) = xx. For B = ( ) as in (3.1 ) and x = u + vi + zj + wij G B, we have 



x = u — (v i + zj + wij) and nrd(a;) = u 2 — av 2 — bz 2 + abw 2 . 

A TLp- lattice of B is a finitely generated Z^-submodule I of B such that FI = B. An 
order O of B is a Zp-lattice which is also a subring of B. A maximal order of B is an order 
which is not properly contained in any other order. Let Oq(1) C B be a maximal order in 
B. 

A right fractional O-ideal is a Zjr- lattice J such that its right order Or(I) = {x G B : 
x/ C J} is equal to 0; left ideals are defined analogously. 

Let K D F be a field containing F. Then S^- = 5 (g)^ is a quaternion algebra over K, 
and we say if splits B if B K = M 2 (K). 

Let f be a noncomplex place of F, and let F v denote the completion of F at v. Then 
there is a unique quaternion algebra over F v which is a division ring, up to isomorphism. 
We say B is unramified (or split) at t> if F v splits -B, otherwise we say B is ramified at t> . 
The set 5 of ramified places of B is a finite set of even cardinality which characterizes B up 
to isomorphism, and conversely given any such set S there is a quaternion algebra over B 
ramified exactly at the places in S. We define the discriminant D of B to be the ideal of 
Tip given by the product of all finite ramified places of B. 



Let 91 C 7j F be an ideal which is coprime to the discriminant D. Then there is an 
isomorphism 

un : O (l) ^ O (l) ®z F = M 2 (Z F , m ) 
where denotes the completion of Zp at 91. Let 

O (91) = {x G Co(l) : is upper triangular modulo 91} ; 

the order C? (9t) is called an Eichler order of level 91. We will abbreviate (9 = (9 (9t). 

We number the real places i>i, . . . , v n of F so that B is split at v i, . . . , v r and ramified at 
v r+ i, . . . ,v n , so that 

B(g) Q R = M 2 (R) r x H n " r 

where EI = ^ — ~|T~~^ * s ^ e division ring of Hamiltonians. If B is ramified at all real places 

(i.e. r = 0) then we say that B is (totally) definite, and otherwise we say B is indefinite. 
The arithmetic properties of the algebra B and its forms are quite different according as 
B is definite or indefinite, and so we consider these two cases separately. Using an adelic 
language, one can treat them more uniformly (though to some extent this merely repackages 
the difference) — we refer to §8 for this approach. 

First, suppose that B is indefinite, so that r > 0. The case B = M 2 (Q) corresponds to 
the classical case of elliptic modular forms; this was treated in §1, so we assume B ^ M 2 (Q). 
Let 

i^-.B^ M 2 {R) r 

denote the map corresponding to the split embeddings v±, . . . , v r . Then the group 

B* = {7 G B x : det^ = (nrd7); > for % = 1, . . . , r} 
acts on W by coordinatewise linear fractional transformations. Let 

Ol = O x HB*. 

Under the assumption that F has strict class number 1, which we maintain, we have 

o\ = z$o? 

where Of = {7 G O : nrd(7) = 1}. Let 

T = r 3 (9I) = l^OI) C GL+(M) r . 

and let Pr = r/Z£. Then Pr is a discrete subgroup of PGL^ (IR) r which can be identified 
with 

pr = U^i x )/{±iKPSL 2 (R). 

Definition 3.2. Let B be indefinite. A quaternionic modular form for B of parallel weight 
2 and level 91 is a holomorphic function f : W — > C such that 

, Q x ,/ \ , faiZi + h a r z r + b r \ ( yr {^Zi + di) 2 \ 

3.3 f^z) = f — , . . . , — =11 77 f(z) 

\c 1 z 1 + d 1 c r z r + d r J det7j J 

for all 7 G r 3 (9l). 
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Analogous to (2.4), we define 



(3-4) (/|7)(*) = /(7*)II 



_det7i_ 



for / : W ->• C and 7 G B£; then ([3^ is equivalent to (/ | 7) (2) = /(*) for all 7 G r^(7t). 

We denote by Mf (7t) the space of quaternionic modular forms for B of parallel weight 2 
and level 7t, a finite-dimensional C-vector space. 

A quaternionic modular form for B = M^-F) is exactly a Hilbert modular form over F; 
our presentation in these three sections has been consciously redundant so as to emphasize 
this similarity. (We could recover the definition of cusp forms given in Section 1 if we also 
impose the condition that the form vanish at the cusps.) As we will see later, this similarity 
is less apparent when the general and more technical theory is exposited. 



The Hecke operators are defined on M^(7t) following their definition in (2.7). Let p be a 
prime of with p { 91, and let p be a totally positive generator of p. Define 

0(p) = 0*\ {vr G + : nrd(7r)Z F = p} = 0*\ {n G + : nrd(vr) = p] . 

The set 0(p) has cardinality iVp + 1. The Hecke operator T p is then given by 

(3-5) (T P f)(z)= Yl (/MO*)- 

7ree(p) 

The set 0(p) admits an explicit description as follows. As above, let F p = Z^/p be the 
residue field of p, and let l p : O <^-> M 2 (Zi? p ) be a splitting. Then the set 0(p) is in bijection 
with the set of left ideals of O by 7r h- > Ott. This set of left ideals is in bijection [HJ Lemma 
6.2] with the set P 1 (F p ): explicitly, given the splitting i p , the left ideal corresponding to 
a={x:y)e P x (F p ) is 

(3.6) J a = 0^{l o) + ° p - 

By strong approximation Theoreme III. 4. 3], each of the ideals J a is principal, so J a = 
07i a with nrd(vr a ) = p for all a G P^Fp). Therefore, we have 6(p) = {n a : a G P^Fp)}. 
This definition reduces to the one given in ( |2.7 ) for Hilbert modular forms with the choices 

VToo = and 7r a = for a G F p . 

Having treated Hilbert modular forms in the previous section, now suppose that B ^ 
M.2(F), or equivalently that B is a division ring. Then a modular form is vacuously a 
cusp form as there are no cusps! We then refer to quaternionic modular forms equally 
well as quaternionic cusp forms and let S-f (9T) = Mf (7t). Here, a cusp form / gives a 
holomorphic differential r-form (27ri) r /(z 1 , . . . , z r ) dz\ ■ ■ ■ dz r on the associated quaternionic 
Shimura variety Xq (71) = r^(7t)\'H T ', a complex variety of dimension r. 

The important case for us will be when r = 1. Then T^(7t) C PGL^(IR) acts on the upper 
half-plane and the quotient r^(7t)\"H can be given the structure of a Riemann surface, known 
as a Shimura curve. In this simple case, a cusp form for B is simply a holomorphic map 

/ : n C such that f(jz) = (c 1 z+d 1 ) 2 f(z) for all 7 G r^(7t), where 71 = ^(7) = r£ 

and t>i is the unique split real place of F . 
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Next, suppose that B is definite, so that r = 0. Recall that O = Oofil) C 0o(l) is 
an Eichler order of level 91. A right fractional (9-ideal is invertible if there exists a left 
fractional CMdeal J -1 such that = O, or equivalently if I is locally principal, i.e., for 
each (finite) prime ideal p of TLp, the ideal I p is a principal right O p -ideal. If / is invertible, 
then necessarily J -1 = {i e B : iJ C £)}. 

Let /, J be invertible right fractional CMdeals. We say that I and J are in the same right 
ideal class (or are isomorphic) if there exists aniG5 x such that / = xJ, or equivalently if 
J and J are isomorphic as right (9-modules. We write [I] for the equivalence class of I under 
this relation and denote the set of invertible right CMdeal classes by CI O. The set CI O is 
finite and H = # CI O is independent of the choice of Eichler order O = OofiX) of level 9T. 

Definition 3.7. Let B be definite. A quaternionic modular form for B of parallel weight 2 
and level 91 is a map 

f : C10 O (9T) -> C. 

The space of quaternionic modular forms M 2 B ( < yt) is obviously a C- vector space of dimen- 
sion equal to H. 

A modular form for B which is orthogonal to the (1-dimensional) subspace of constant 
functions is called a cusp form for B; the space of such forms is denoted S% (OT). 

The Hecke operators are defined on M 2 B ( i yi) as follows. Let p be a prime ideal of with 
p | 9T. For a right (9-ideal I with nrd(J) coprime to p, the Hecke operator T p is given by 

(3-8) (T P f)([I])= E W])' 

JCJ 
nrd(J/- 1 )=p 



the sum over all invertible right CMdeals JCI such that nrd(J) = pnrd(J). As in (3.5), 
this sum is also naturally over P x (Fp), indexing the ideals of norm index p. This definition 
does not depend on the choice of representative I in its ideal class and extends by linearity 
to all of Si {VI). 

Consequent to the definitions in the previous paragraphs, we may now consider the Hecke 
modules of quaternionic cusp forms over F for the different quaternion algebras B over 
F. These spaces are related to each other, and thus to spaces of Hilbert modular forms, 
according to their arithmetic invariants by the Jacquet-Langlands correspondence as follows. 

Theorem 3.9 (Eichler-Shimizu- Jacquet-Langlands). Let B be a quaternion algebra over F 
of discriminant T) and let be an ideal coprime to D . Then there is an injective map of 
Hecke modules 

s?{m) s 2 {T>m) 

whose image consists of those Hilbert cusp forms which are new at all primes p | D. 

Proof. See Jacquet and Langlands [36, Chap. XVI], Gelbart and Jacquet [271 §8] an d work 
of Hida [32J; another useful reference is Hida [331 Proposition 2.12] who deduces Theorem 



3.9 from the representation theoretic results of Jacquet and Langlands. □ 



Theorem 3J3 yields an isomorphism 

si(m) s 2 (m) 

when the quaternion algebra B has discriminant D = (1). Since a quaternion algebra must 
be ramified at an even number of places, when n = [F : Q] is even we can achieve this for 
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the definite quaternion algebra B which is ramified at exactly the real places of F (and no 
finite place). When n is odd, the simplest choice is to instead take B to be ramified at all 
but one real place of F (and still no finite place), and hence B is indefinite (and g = 1). 

Remark 3.10. Note that in general a space of newforms can be realized as a Hecke mod- 
ule inside many different spaces of quaternionic cusp forms. Indeed, for any factorization 
0Jt = Qyi with D squarefree and coprime to D, we consider a quaternion algebra B of 
discriminant D (ramified at either all or all but one real place of F) and thereby realize 
Sif (9?) = 5 , 2(9Jl) s,_new . For example, if p, q are (rational) primes, then the space S , 2(pg) new 
of classical newforms can, after splitting off old subspaces, be computed using an indefinite 
quaternion algebra of discriminant 1 or pq (corresponding to a modular curve or a Shimura 
curve, respectively) or a definite quaternion algebra of discriminant p or q. 

Our main conclusion from this section is that to compute spaces of Hilbert cusp forms it 
suffices to compute instead spaces of quaternionic cusp forms. The explicit description of 
S'.f (?Tt) as a Hecke module varies according as if B is definite or indefinite. 



4. Definite method 

In this section, we discuss a method for computing Hilbert modular forms using a definite 
quaternion algebra B. We continue with our notation and our assumption that F has strict 
class number 1. We accordingly call the method in this section the definite method: it is 
a generalization of the method of Brandt matrices mentioned briefly in §2 and was first 
exhibited by Eichler [22] and Pizer jlH] for F — Q, but the first explicit algorithm was given 
by Socrates and Whitehouse [60] . 

Let Ji, . . . , Ih be a set of representative right ideals for CI O, with H = # CI O. As vector 
spaces, we have simply that 

Mf (9t) = Map(ClC,C) ®f =1 CIi, 

associating to each ideal (class) its characteristic function. Let 0\ = Ol{Ii) be the left order 
of and let a = #(0*/Z*). 



The action of the Hecke operators is defined by (3.8): we define the pth-Brandt matrix 



T(p) for O to be the matrix whose (i,j)th entry is equal to 

(4.1) b(p) id = #{J C Ij : nrd(J/7 1 ) = p and [J] = [/*]} G Z. 

The Brandt matrix T(p) is an H x if-matrix with integral entries such that the sum of the 
entries in each column is equal to iVp + 1. The Hecke operator T p then acts by T(p) on 
; C/j (on the right), identifying an ideal class with its characteristic function. 

The Brandt matrix is just a compact way of writing down the adjacency matrix of the 
graph with vertices X = CI O where there is a directed edge from Ij to each ideal class which 
represents an ideal of index iVp in Ij. Indeed, consider the graph whose vertices are right 
(9-ideals of norm a power of p and draw a directed edge from I to J if nrd(J7 -1 ) = p. Then 
this graph is a k- regular tree with iVp + 1 edges leaving each vertex. The above adjacency 
matrix is obtained by taking the quotient of this graph by identifying two ideals if they are 
in the same ideal class. 

Alternatively, we may give an expression for the Brandt matrices in terms of elements 
instead of ideals. A containment J C Ij of right CMdeals with [Ij] = [J] corresponds to 
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an element x G JIf C Ijlf 1 via J = xli, and we have nrd (JJ^ 1 ) = p if and only if 
nrd(x)Zi? = p. 

Writing JI^ 1 = xOi, we see that x is unique up to multiplication on the right by Of. 
We have Of = (O^Z* and (O^ n Z* = {±1}, so 2e { = #(O i ) 1 x . To eliminate the 
contribution of the factor Z£, we normalize as follows: let p be a totally positive generator 
for p and similarly for nrd(ij) for i = 1, . . . , H. Then x G ij/i -1 can be chosen so that 
md(x)(qj/qi) = p and is unique up to multiplication by . 

(4.2) 6(p) M = -L# {* G J,/- 1 : nid(s)| = . 



The advantage of the expression (4.2) is that it can be expressed simply in terms of a 
quadratic form. Since B is definite, the space B c — >■ B ®q R = HP = IR 4 ™ comes equipped 
with the positive definite quadratic form Trnrd : B — > R. If J is a Z^-lattice, then J = Z 4n 
embeds as a Euclidean lattice J e — >■ lR 4n with respect to this quadratic form. It follows that 
one can compute 6(p)ij by computing all elements x G -fj-^ -1 such that Tr^/g^) nrd(x) < 
Trp, a finite set. 

Before giving references for the technical details about how the Brandt matrices above are 
computed explicitly, we pause to give three examples. 

f — \ — 23 \ 

Example 4.3. Consider the quaternion algebra B = I — '— J, so that B is generated by 



i,j subject to i 2 = —1, j 2 = —23 and ji = —ij. We have the maximal order O = Oq(1) = 
7h@ r LiQ>7Lk@ r Lik where k = (j + l)/2. We consider the prime p = 2; we have an embedding 



O ^ M 2 (Z 2 ) where i, k ^ ^ Jj , °J (mod 2). 



We begin by computing the ideal classes in O. We start with G\ = [O]. We have 3 ideals 
of norm 2, namely J (0:1) = 20 + ikO, = 20 + {i + l)kO, and J (lj0 ) = 20 + jfeO. If one of 
these ideals is principal, then it is generated by an element of reduced norm 2. The reduced 
norm gives a quadratic form 

nrd : O -> Z 

x + + + wik \— > x 2 + xz + y 2 + yu> + 6^ 2 + 6u> 2 

We see immediately that nrd(x + yi + zk + wik) = 2 if and only if z = w = and x = y = ±1. 
This shows that I(ia) = (i + 1)0 is principal but J(i :0 ) and J( 0: i) are not. Note also that we 
find 2ei = 4 solutions matching (4.2). We notice, however, that z/(i : o) = J(o:i)> so we have 
just a second ideal class C 2 = [/(o-i)J- 

Now of the two ideals contained in /( 0;1 ) of norm 4, we have iJJ^ = AO + i(k + 2)0 belongs 
to Ci whereas = ^ + + 2k + ik)0 gives rise to a new ideal class C 2 = [ipa)]- If 

we continue, each of the two ideals contained in -^o-i) °f norm 8 belong to C±, and it follows 
that H = #C1£> = 3. 

/I 1 N 

From this computation, we have also computed the Brandt matrix T(2) =12 13 

V° 1 °y 

Indeed, the first column encodes the fact that of the three right CMdeals of reduced norm 
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2, there is one which is principal and hence belongs to C\ and two that belong to Ci- We 
think of this matrix as acting on the right on row vectors. 

The characteristic polynomial of T(2) factors as (x — 3)(x 2 + x — 1). The vector (1, 1, 1) is 
an eigenvector with eigenvalue 3 which generates the space of constant functions and gives 
rise to the Eisenstein series having eigenvalues a p = p + 1 for all primes p ^ 23. The space 
of cusp forms S 1 ^ (1) is correspondingly of dimension 2 and is irreducible as a Hecke module. 
The Hecke module 5^(1) can be explicitly identified with S l 2(23) using theta correspondence. 
For example, the series 

Q^q) = £ <f rd « = g - 2 +-+f 2 +^+6, W = 1+4q 2 + 4g 4 + 4g8 + gg 10 + 

is the g-expansion of a modular form of level 23 and weight 2 and corresponds to (the 
characteristic of) C\. For more details, we refer to Pizer [48, Theorem 2.29], where the first 
computer algorithm for computing Brandt matrices over Q is also described. 

Now we give an example over a quadratic field. 

Example 4.4. Let F = Q(\/5). Then Z F = Z[w] where w = (l + v / 5)/2 satisfying w 2 — w — 1 = 
0. The prime 61 splits in F; let 9T = (3w + 7)Z^ be one of the primes above it. We consider 

the (Hamilton) quaternion algebra B = ( — j- — j over F of discriminant D = (1). We 



have the maximal order 

O (l) =Z F ®iZ F ®kZ F ®ikZ F , where k = ± '— 

and the Eichler order O C C?o(l) of level 91 given by 

C = Z F © (3w + 7)iZ F © (— 30i + A;)Z F © [w + 20z + zfc)Z F . 

The class number of O is H = 3. The following ideals give a set of representatives for 
CI O: we take h = 0, 

I 2 = 20 + ((w + 2) - (2w + 2)i + (-1 + 3w)ik)0 

= 2Z F © (6w + U)iZ F © ((w + 1) + (-2«; + 5)i - k)Z F © (1 - 6i + wiA;)Z F 

and J 3 = 20 + ((to + 1) + (1 - w)i + (2 - 2w)A;)0. 

We compute the orders = #(Oi/Z F ) as ex = 2, e 2 = 5 and e 3 = 3. For example, the 
element u = {2w)i — k — wik G O satisfies the equation u 2 + 1 = 0, and so yields an element 
of order 2 in O x fL F . It follows that none of these orders are isomorphic (i.e., conjugate) in 
B. 

The first few Brandt matrices are: 

(1 5 3\ 

T(2) = 2 , T(v / 5) = 
\2 2/ 




T(w + 3) = 4 2 I , T(w-A) 
\4 6 

We note that iV(w + 3) = A^(w - 4) = 11. 
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For example, the first column of the matrix T(2) records the fact that of the 5 = N(2) + 1 
right CMdeals of norm (2), there is exactly one which is principal, two are isomorphic to I2 
and the other two are isomorphic to 1%. 

The space Sj (OT) of cusp forms is an irreducible 2-dimensional Hecke module, represented 
by a constituent form / with corresponding eigenvector (2, 5u>, — 3w — 3). In particular, the 
ring of Hecke operators T/ = Z[a p ] restricted to / is equal to 1/ = Z[w], by coincidence. We 
have the following table of eigenvalues for /. 



p 


(2) 


iw + 2) 


(3) (w + 3) (w-A) 


iVp 


4 


5 


9 11 11 


<h(f) 


2w - 2 


-3w + 1 


— w — 2 4u> — 2 — w 



For further discussion of the geometric objects which arise from this computation, see the 
discussion in Section [6] 

Finally, one interesting example. 



Example 4.5. Let F = Q(v 15) and let 9T = (5, y 15). Then there exists a cusp form con- 
stituent of an irreducible space of dimension 8 in S^OT) such that no single Hecke eigenvalue 
generates the entire field E of Hecke eigenvalues. 



E 



gal 



17, i, A y/u) 
17, i, v^«) 




Here, w = (5 + vl7)/2 and i£ ga i is the Galois closure of -E. Each Hecke eigenvalue a p for 
this form / generates a proper subfields of E. (There are also examples of this phenomenon 
over Q, and they are related to the phenomenon of inner twists; this was analyzed over Q 
by Koo, Stein, and Wiese |4~4"].) 

With these examples in hand, we now give an overview of how these computations are 
performed; for more detail, see work of the first author [15]. It is clear we need several 
algorithms to compute Brandt matrices. First, we need a basic suite of algorithms for 
working with quaternion orders and ideals; these are discussed in work of Kirschmer and the 
second author jH] Section 1], and build on basic tools for number rings by Cohen [7j. As 
part of this suite, we need a method to compute a maximal order, which is covered by work 
of the second author [70]. Next, we need to compute a set of representatives for GIO and 
to test if two right (9-ideals are isomorphic: this is covered by Kirschmer and the second 
author [HI Sections 6-7] , including a runtime analysis. To compute a set of representatives, 



we use direct enumeration in the tree as in Example 4.3 and a mass formula due to Eichler 
as an additional termination criterion. To test for isomorphism, we use lattice methods to 
find short vectors with respect to the quadratic form Tr nrd. 
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In this method, to compute with level 9T one must compute a set of representatives CI O = 
CI Oq anew. The first author has given an improvement on this basic algorithm, allowing 
us to work always with ideal classes Cl(9o(l) belonging to the maximal order at the small 
price of a more complicated description of the Hecke module. The proof of correctness for 
this method is best explained in the adelic language, so we refer to Section 8 for more detail. 

Let Ix,...,Ih be representatives for Cl(9o(l) and suppose that 9? is relatively prime to 
nrd(/j) for each i — this is made possible by strong approximation. As before, let p be a prime 
of Z F with p \ D. Let p be a totally positive generator for p, and let qi be a totally positive 
generator for nrd(Zj). By our notation, we have Oi(Ii) = Oo(l);- Let Tj = O (1)* jZ F . For 
each i, j, consider the set 



0(p) M = r,-\ j* E Ijlr 1 : nrd(x)^ = pj 



where Tj acts by the identification O (l)f jZ F = O (l)i / {±1}- Via a splitting isomorphism 

Uft : O (l) ^ Oo(l) ® Zjr.gt = M 2 (Z Fi ^), 

the group 0o(l) x acts on P^Z^/D"?) and since 0o(l) ® Z^ = Oo(l)j ® Zi?^ for each % (since 
nrd(Ji) is prime to 9t), the group F t = O {l)* /Z F similarly acts on F 1 (Z F /m). 
We then define a Hecke module structure on © i=1 C[r i \P 1 (Z f '/9 r l)] via the map 

ciTj-XP^Zjr/m)] c[r i \p 1 (z J? /m)] 

on each component. It is a nontrivial but nevertheless routine calculation that this Hecke 
module is isomorphic to the Hecke module M 2 defined by the Brandt matrices at the 
beginning of this section. 



Example 4.6. We keep the notations of Example 4.4 The quaternion algebra B has class 
number 1, thus the maximal order O (l) is unique up to conjugation. We have T = 
Oq(1) x /Zp has cardinality 60. We consider the splitting map 

l m : O (l) O (l) ®z F (Z F /(H) = M 2 (Z F /9t) 

given by 

_ (11 0\ _ /47 58\ _ r . /29 28\ 

^ W = 1,37 50 J ' L ^ k) = \l% 33 ) ' = 1,16 14 ) " 

(One directly verifies that %(ifc) = We let T act on ¥ 1 {Z F / ( yC) on the left via 

I<n- By the above discussion, we have 

Mf(at) ^ cpvp^Zf/^)]. 

The action of T on P 1 (Z^/9 f t) has three orbits which are represented by x\ = (1 : 0), 
x 2 = (1 : 1) and x 3 = (23 : 1) whose stabilizers have cardinality t\ = 2, e 2 = 5 and e 3 = 3. 
Thus M^i^yi) is a free module generated by the orbits Txi, Tx 2 and Tx 3 . Writing down the 



Hecke action in that basis, we obtain the same Hecke operators as in Example 4.4 



Remark 4.7. The approach presented above has some advantages over the usual definition 
of Brandt matrices above. First of all, it is better suited for working with more general 
level structures, such as those that do not come from Eichler orders. For example, adding a 
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character in this context is quite transparent. Secondly, when working over the same number 
field, a substantial amount of the required data can be precomputed and reused as the level 
varies, and consequently one gains significantly in the efficiency of the computation. 



5. Indefinite method 

In this section, we discuss a method for computing Hilbert modular forms using a indefinite 
quaternion algebra B with r — 1. We accordingly call our method the indefinite method. 
The method is due to Greenberg and the second author [30]. We seek to generalize the 
method of modular symbols by working with (co)homology. We continue to suppose that F 
has strict class number 1, and we assume that B ^ M 2 (Q) for uniformity of presentation. 

Recall that in this case we have defined a group Y = r^(9T) C GLj(R) such that 
Pr = T/Zp C PGI4(M) is discrete; the quotient X = X^(m) = T\H is a Shimura curve 
and quaternionic cusp forms on B correspond to holomorphic differential 1-forms on X. Inte- 
gration gives a Hecke-equivariant isomorphism which is the analogue of the Eichler-Shimura 
theorem, namely 

S*(<R) © SfW) ^ H\X$ (91), C). 

We recover Sjf (Ot) by taking the +- eigenspace for complex conjugation on both sides. Putting 
this together with the Jacquet-Langlands correspondence, we have 

s 2 {®mf~ ncw = s^(m) = H\x^{m),c) + - 

We have the identifications 

H\X, C) = H\X^{m), C) S H\rfl{Vl), C) = Hom(r£(9t), C) = Hom(r, C). 
To complete this description, we must relate the action of the Hecke operators. Let p \ 



be prime and let p be a totally positive generator of p. As in (3.5), let 
(5.1) 0(p) = 0*\ {tt G 0+ : nrd(Tr) = p} 

and choose representatives ir a for these orbits labeled by a G P X (F P ). Then any 7 G V by 
right multiplication permutes the elements of 0(p), and hence there is a unique permutation 
7* of P X (F P ) such that for all a G P X (F P ) we have 

7Ta7 = 5a7T 7 *a 

with 5 a G T. For / G H 1 ^, C) = Hom(r, C), we then define 

(5-2) (V)( 7 )= £ f(Sa)- 

7r a eG(p) 

In a similar way, we compute the action of complex conjugation via the relation 

(T 00 /)( 7 ) = /(5) 

where 7x7 = 5/j, and \i G O x \ 0\ . 

We begin with two examples, to illustrate the objects and methods involved. 



Example 5.3. Let F = Q(y29). Then Z,p = Z,[w] is the ring of integers of F where w = 
(1 + y / 29)/2 satisfies w 2 — w — 7 = 0. Indeed F has strict class number 1 and u = w + 2 is 
a fundamental unit of F. 
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Let B = ( ^ J , so B is generated over F by i,j subject to i 2 = —1, j 2 = u, and 

ji = —ij. The algebra B is ramified at the prime ideal 2Zp and the nonidentity real place 
of B (taking \/29 — \/29) and no other place. The identity real place gives an embedding 



ioo : B e — >■ B <^>p m. = ivi2 ^ 

A / \/u \ Z2.27... 



z ' J ^ l -f y ' V V -2.27.. 



Let 



O = O (l) = 7L F ® Z F i © Z F j © Z F k 



where k = (l + i)(w + l+j)/2. Then O is a maximal order of B with discriminant D = 27L F . 
Let T = LooiO*) C SL 2 (M) be as above and let X = X(l)c = r\"H be the associated Shimura 
curve. 

Although they are not an intrinsic part of our algorithm, we mention that the area of X 
(normalized so that an ideal triangle has area 1/2) is given by 



A = (^4 /2 C,(2)*(f) = p^fr(2)(4 - 1) = \ 



where $(£•) = n p |2)(^P + 1)> anc i ^ ne genus of X is given by the Riemann-Hurwitz formula 

as 



where e q is the number of elliptic cycles of order q G Z> 2 in V. An explicit formula for e q 
given in terms of class numbers and Legendre symbols yields e 2 = 3 and e q = for q > 2. 
Thus 2<7 — 2 = 0, so (7 = 1. For more details on these formulas and further introduction, see 
work of the second author |68J and the references given there. 

Next, we compute a fundamental domain for T, yielding a presentation for T; we consider 
this as a black box for now. The domain, displayed in the unit disc, is as follows. 
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We obtain the presentation 

T = (7, 7\ *i, 62, 5 3 : 81 = 5 2 2 = 51 = [ 7 , 7IW3 = 1) 

where 

7 = — w — 1 + i — 2j + k 
7' = 2 + 2z + (w - l)j -(w- l)k 
5 X = (2w + 2) +wi + j + 4k 

5 2 = i 

5 3 = (w + 1) + (2w + 3)i -j-k. 
The above method gives the isomorphisms of Hecke modules 

S 2 {2Z F y™ = S*{1) = {f:H^C: f{gz) d{gz) = f(z) dz for all g E T} = H\T, C) + . 

We compute that if 1 (r,C) + = Hom(r,C) + = Cf where / is the characteristic function of 
7, i.e., /( T ) = 1, /(f) = and f(5 t ) = for i = 1, 2, 3. 

We compute the Hecke operator T p for p odd according the definition (5.2). Let p = 
(w + 1)Z F . Then iVp = p — 5. We compute the action of T p on H l (T, C) + given by 
T p f = a p (f)f. The Hecke operators act as a sum over p + 1 left ideals of reduced norm p, 
indexed by F 1 (¥ 5 ). Let 

t p : O ^ M 2 (Z FiP ) = M 2 (Z 5 ) 

Specifying the images modulo p gives them uniquely, as they lift to M 2 (Z 5 ) using Hensel's 
lemma. (Note that j 2 = u = 1 (mod w + 1).) 
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Let Jqo, Jq, . . . , J a be defined by 



Ja J(x:y) 



as in (3.6). Then J a = On a are principal left CMdeals by strong approximation. For example, 
J = Oii - ij) + (w + 2)0 = O-kq where tt = (-it; + 3) + wi + j + ij. 

We compute the Hecke operators as in (5.2). For / : T — » C and 7 G T, we compute 



elements <5 a G T indexed by a G P 1 (F p ) and 7* a permutation of P 1 (F p ) such that 7r a 7 
for all a G P 1 (F p ); then 

(V)(7)= £ /(*«)■ 

The contribution to the sum for / simply counts the number of occurrences of 7 in the 
product 7r a 77r 7 * a = 5 a G T. Carrying out this computation for various primes, we obtain the 
following table. 



iVp 


5 


7 


9 


13 


23 


29 


53 


59 


67 


71 


83 




1 


-2 


5 


-1 


-6 


10 


-1 


10 


8 


-8 


13 


Np + 1 - a(p) 


5 


10 


5 


15 


30 


20 


55 


50 


60 


80 


70 



Here we list only the norm of the prime as the eigenvalue does not depend on the choice 
of prime p of the given norm. This suggests that / corresponds to a base change of a form 
from Q. 

So we look through tables of elliptic curves over Q whose conductor is divisible only by 2 
and 29. We find the curve E labelled 1682cl, where 1682 = 2 • 29 2 , given by 



E:y 2 



xy 



x 3 + x 2 



51318x - 4555676. 



Let Ep denote the base change of E to F. We compute that the twist E' F of E F by — u 
given by 



29, 



x 3 + (-117 + l)x 2 + (-llw - 20)x + (23w + 52) 



E' F : y 2 + (w + l)xy 

has conductor 2Zj7. Since the extension F/Q is abelian, by base change theorems we know 
that there exists a Hilbert cusp form associated to E' F over F which is new of level 2 r Lp, 
which therefore must be equal to /. This verifies that the Jacobian J(l) of X(l) is isogenous 
to E' F . We verify that #E' F (W p ) = Np + 1 — a(p) and (as suggested by the table) that E has 
a 5-torsion point, (—1, — 2w — 5) (and consequently so too does </(!)). 



Example 5.4. To illust rate the Jacquet-Langlands correspondence (Theorem 3.9) in action 
we return to Example 



4.4 



w. 



Recall F = Q(a/5) and w — (1+ yfb)/2. We find the quaternion 
( 3 m I 7*1 *\ 

which is ramified at = (3w + 7)Zjr, a prime of norm 61, 



algebra B = 

V r J 

and one infinite place. The order O = Oq(1) with Z^-basis k,ik, where k = ((w + 1) 
wi+j)/2, is maximal. As above, we compute that A(X) = 10 and g(X) = 2 = dim 5*^(1) 
dim ^(OT) new . Now we have the following fundamental domain: 
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We obtain the presentation 

r(l) = (71, in 72, 72 : [7i, 7i] [72,72] = !)■ 
(In particular, T(l) is a torsion-free group.) We have dimH l (X, C) = 4, and on the basis 
of characteristic functions given by 7i,7i,72,72> the action of complex conjugation is given 



by the matrix 



(1 


1 


\ 





-1 








1 


-1 




-1 


-1 0) 


+ = 


Si(l) 


. Comput 



: this is computed in a way We thus obtain a basis for 



as in Example 4.4 
H\T 2 = 



that Sn(l) is irreducible as a Hecke module, and find that 



H I T, 



w+3 



H 




6 


') 


-4 


-6 



-5 - 
3 4 

if I T„ 



^ I 7!s 



1 



Happily, the characteristic polynomials of these operators agree with those computed using 
the definite method. 

We now give an overview of how these computations are performed: for more details, see 
the reference by Greenberg and the second author [30]. To compute effectively the systems 
of Hecke eigenvalues in the cohomology of a Shimura curve, we need several algorithms. 
First, we need to compute an explicit finite presentation of T with a solution to the word 
problem in T, i.e., given 5 G T, write S as an explicit word in the generators for T. Secondly, 
we need to compute a generator (with totally positive reduced norm) of a left ideal I C O. 
The first of these problems is solved by computing a Dirichlet domain; the second is solved 
using lattice methods. We discuss each of these in turn. 
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Let p G H have trivial stabilizer r p = {g G T : gp = p} = {1}. The Dirichlet domain 
centered at p is the set 



D(p) = {zeH: d(z,p) < d(gz,p) for all g G T} 



where d denotes the hyperbolic distance. In other words, we pick in every T-orbit the closest 
points z to p. The set D(p) is a closed, connected, hyperbolically convex fundamental domain 
for T whose boundary consists of finitely many geodesic segments, called sides, and comes 
equipped with a side pairing, a partition of the set of sides into pairs s, s* with s* = g(s) for 
some g G T. (We must respect a convention when a side s is fixed by an element of order 2, 
considering s to be the union of two sides meeting at the fixed point of g.) 

The second author has proven [S7j that there exists an algorithm which computes a Dirich- 
let domain D for T, a side pairing for D, and a finite presentation for T with a minimal set 
of generators together with a solution to the word problem for T. This algorithm computes 
D inside the unit disc T>, and we consider now T acting on D by a conformal transforma- 
tion which maps p — y 0. We find then that T> can be computed as an exterior domain, 
namely, the intersection of the exteriors of isometric circles 1(g) for elements g G T, where 

1(g) = {z G C : \cz + d\ = 1} where g = ^ G SU(1, 1). We find generators for the 

group by enumerating short vectors with respect to a certain positive definite quadratic form 
depending on p. 

The solution to the word problem comes from a reduction algorithm that finds elements 
which contribute to the fundamental domain and removes redundant elements. Given a 
finite subset G C V \ {1} we say that 7 is G-reduced if for all g G G, we have <i(70,0) < 
d(g , y0,0). We have a straightforward algorithm to obtain a G-reduced element, which we 
denote 7 H- redc(7): if ^(70, 0) > d(gj0, 0) for some g G G, set 7 := #7 and repeat. When 
the exterior domain of G is a fundamental domain, we have rede (7) = 1 if and only if 7 G I\ 
This reduction is analogous to the generalized division algorithm in a polynomial ring over 
a field. 

Example 5.5. Let F be the (totally real) cubic subfield of Q(Ci3) with discriminant dp = 169. 
We have F = Q(b) where 6 3 + 4b 2 + 6 — 1 = 0. (F has strict class number 1.) 

The quaternion algebra B = ( — — j has discriminant D = (1) and is ramified at 2 of 



F 

the 3 real places of F. We take O to be an Eichler order of level p = (b + 2), a prime ideal 
of norm 5; explicitly, we have 

C = /.^p © (0 + 2)2Zp © L F © iLp. 



We compute a fundamental domain for the group T = Vq (p) = ^(Ol) / {±.1} . We take 

p = 9/lOz en. 

We first enumerate elements of O by their absolute reduced norm. Of the first 260 ele- 
ments, we find 29 elements of reduced norm 1. 
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Let G be the set of elements which contribute to the boundary. For each g ^ G, we 
compute redct^)- Each in fact reduces to 1, so we are left with 8 elements. 
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We next enumerate elements in O moving the center in the direction of the infinite vertex 
vj. We find an enveloper, an element g such that v-/ lies in the interior of 1(g), and reduce 
to obtain the following. 




The domain now has finite area. We now attempt to pair each vertex to construct a side 
pairing. 




For example, the first vertex v\ pairs with v$, pairing the isometric circle I(gi) of g± with 
I(9t)- We continue, but find that Vg does not pair with another vertex. 
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Indeed, gg(v 9 ) does not lie in the exterior domain. We compute the reduction red^^g; v 9 ), 
analogously defined. We then obtain a domain with the right area, so we are done! 




(In order to get an honest side pairing, we agree to the convention that a side 1(g) which 
is fixed by an element g, necessarily of order 2, is in fact the union of two sides which meet 
at the unique fixed point of g; the corresponding vertex then appears along such an side.) 

As to the second problem, that of principalizing ideals, we refer to work of Kirschmer and 
the second author [4TJ §4], which relies again upon enumeration in a lattice. This method can 
be improved using the computation of a fundamental domain, which allows one to further 
effectively reduce lattice elements. 

As in the previous section, when the level varies, by changing the coefficient module 
one can work always with the group T(l) associated to the maximal order Oq(1) rather than 
recomputing a fundamental domain for T = ro(9T) each time. This is a simple application 
of Shapiro's lemma: the isomorphism 

H\r, c) tf^rxi), Coind£ (1) c), 
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is an isomorphism of Hecke modules: we give Coindp V ; C = Hom(r(l)/r,C) the natural 
structure as a r(l)-module, and for a cochain / G H 1 ^, Coindp C), we define 

(5-6) (V)(7)= £ W 

as in (5.2). We require still that O(p) C + ; after enumerating a set of representatives for 



T\r(l), we can simply multiply elements ir a G Oq(1) by the representative for its coset. 

To conclude, we compare the above method to the method of modular symbols introduced 
when F = Q. The Shimura curves X = Xg(Vl) do not have cusps, and so the method of 
modular symbols does not generalize directly. However, the side pairing of a Dirichlet domain 
for T gives an explicit characterization of the gluing relations which describe X as a Riemann 
surface, hence one obtains a complete description for the homology group Hi(X,I*). Paths 
are now written {v , jv} for v a vertex on a side paired by 7 G G. The analogue of the Manin 
trick in the context of Shimura curves is played by the solution to the word problem in T, 
and so in some sense this can be seen as a partial extension of the Euclidean algorithm to 
totally real fields. Our point of view is to work dually with cohomology, but computationally 
these are equivalent [301 §6]- 



6. Examples 

In this section, we provide some detailed examples to illustrate applications of the algo- 
rithms introduced in the previous sections. 

Our first motivation comes from the association between Hilbert modular forms and 
abelian varieties. The Eichler-Shimura construction (see Knapp [38, Chap. XII] or Shimura 
[591 Chapter 7]) attaches to a classical newform f(z) = Y^ n a n q n G S2(N) new sua abelian vari- 
ety Af with several properties. First, Af has dimension equal to [Ef : Q] where Ef = Q({a„}) 
is the field of Fourier coefficients of / and Endq(Af) <g> Q = Ef. Second, Af has good reduc- 
tion for all primes p \ N . Last, we have an equality of L-functions 

L(A f ,s)= J] L(r,s) 

creHom(£; / ,C) 

where f a is obtained by letting o act on the Fourier coefficients of /, so that 

n=l 

The abelian variety Af arises as the quotient of Jq(N) by the ideal of End Jq(N) generated 
by T n - a n . 

The simplest example of this construction is the case of a newform / with rational Fourier 
coefficients Q({a n }) = Q. Then A f is an elliptic curve of conductor N obtained analytically 
via the map 

X (N)^A f 

t 1 y z r = 2ni / f(z) dz 

J ioo 
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00 

n 

n=l 



The equality of L-functions is equivalent to the statement #Af(W p ) = p+l — a p for all primes 
p, and this follows in the work of Eichler and Shimura by a comparison of correspondences 
in characteristic p with the Frobenius morphism. 

The Eichler-Shimura construction extends to the setting of totally real fields in the form 
of the following conjecture (see e.g. Darmon [TTj, Section 7.4]). 

Conjecture 6.1 (Eichler-Shimura). Let f G S2(91) ncw be a Hilbert newform of parallel weight 
2 and level 91 over F. Let Ef = Q({a n }) be the field of Fourier coefficients of f . Then there 
exists an abelian variety Af of dimension g = [Ef : Q] defined over F of conductor 91 9 such 
that End F (A f ) <g> Q ^ E f and 

L(A f ,s)= J] L(r,s). 

creHom(_B / ,C) 



In particular, if / has rational Fourier coefficients, Conjecture |6.1| predicts that one can 
associate to / an elliptic curve Af defined over F with conductor 91 such that L(Af,s) = 
L(f,s). 



Conjecture 6.1 is known to be true in many cases. It is known if / appears in the coho- 
mology of a Shimura curve (as in Sections 3 and 5 above; see also Zhang [74]): this holds 
whenever n = [F : Q] is odd or if 91 is exactly divisible by a prime p. If / appears in the 
cohomology of a Shimura curve X and / has rational Fourier coefficients, then when the 
conjecture holds there is a morphism X — Y Af and we say that Af is uniformized by X. 



However, Conjecture |6.1| is not known in complete generality: by the Jacquet-Langlands 
correspondence, one has left exactly those forms / over a field F of even degree n and 
squarefull level (i.e. if p | 91 then p 2 | 91), for example 91 = (1). One expects still that the 
conjecture is true in this case: Blasius pQ, for example, shows that the conjecture is true 
under the hypothesis of the Hodge conjecture. 



The converse to the Eichler-Shimura construction (Conjecture 6.1 ) over Q is known as the 
Shimura- Taniyama conjecture. An abelian variety A over a number field F is of GL2-type 
if Endi? A ® Q is a number field of degree dim A. The Shimura- Taniyama conjecture states 
that given an abelian variety A of GL2-type over Q, there exists an integer N > 1 and a 
surjective morphism Jq(N) — > A. This conjecture is a theorem, a consequence of the proof 
of Serre's conjecture by Khare and Wintenberger In the setting of totally real fields, 
the analogous conjecture is as follows. 

Conjecture 6.2 (Shimura- Taniyama). Let A be an abelian variety of GL2-type over a totally 
real field F . Then there exists a Hilbert newform f of parallel weight 2 such that Ef = 
End F (A) <g> Q and 

L(A,s)= H L(r,s). 

ueHomfEj.C) 

If both of these conjectures are true, then by Faltings' isogeny theorem, the abelian variety 
A in Conjecture 6.2 would be isogenous to the abelian variety Af constructed in Conjecture 



6.1| An abelian variety of GL 2 -type over a totally real field F which satisfies the conclusion 



of Shimura- Taniyama conjecture (Conjecture 6.2) is called modular. 



There has been tremendous progress in adapting the techniques initiated by Wiles [TTj 
which led to the proof of the Shimura- Taniyama conjecture in the case F = Q — to totally real 
fields. For example, there are modularity results which imply that wide classes of abelian 
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varieties of GL2-type are modular: see for example work of Skinner and Wiles [6T| [62] , 
Kisin 1121133], Snowden [63] and Geraghty [28], and the results below. However, a complete 



proof of Conjceture 6.2 remains elusive. 



With these conjectures in mind, we turn to some computational examples. We begin with 
F = Q(v / 5) and let w = (1 + v / 5)/2. For some further discussion on this case, see work of 
the first author [14|. 

Example 6.3. We consider the spaces S^^T) with iV(9T) < 30. 

We compute using the definite method of Section 4 using the quaternion algebra over F 
which is ramified at no finite place. We find that dim ^(91) = for all ideals 9T of with 
N{% < 30 and dimS^Dt) = 1 for 01 with norm 31. 

Now let 91 be either of the primes above 31. By the Jacquet-Langlands correspondence 



(Theorem 3.9) and the accompanying discussion, the space S^OT) can also be computed using 
the indefinite quaternion algebra B ramified at 01 and one real place. We find agreeably that 
the Shimura curve Xq (1) has genus one and thus its Jacobian is an elliptic curve. By a 
search, we find the elliptic curve 

A: y 2 + xy + wy = x 3 — (1 + w)x 2 

of conductor 91 = (5 + 2w) of norm 31. 

We now show that A is modular. In brief, we verify that the Galois representation p A 3 : 
Gal(F/F) -)• GL 2 (F 3 ) C GL 2 (F 9 ) is irreducible, solvable, and ordinary since #A(F 9 ) = 8 
(so 3 \ a 3 (A) = 9 + 1 — 8 = 2), so A is modular by Skinner and Wiles [HI]. Since A has 
conductor p and the space S^p) has dimension 1, it follows by Falting's isogeny theorem 
that A is isogeneous to the Jacobian of the Shimura curve Xq (1). 

In particular, the elliptic curves over F = Q(\/5) with smallest conductor have prime 
conductor dividing 31 and are isogenous to A or its Galois conjugate. 

Example 6.4. Let 91 C Z F be such that iV(9tj < 100. 

If 91 { 61, then dim S2(9T) new < 2: each newform / of conductor 91 has rational Fourier 
coefficients and is either the base change of a classical modular form over Q or associated to 
an elliptic curve which is uniformized by a Shimura curve. 

We compute each space in turn using the definite method. Those forms / which have the 
property that %,(/) = cip(/), where ~~ denotes the Galois involution, are candidates to arise 
from base change over Q. For each of these forms /, we find a candidate classical modular 
form g using the tables of Cremona [9] : such a form will have conductor supported at 5 and 
the primes dividing iV(9T). Each curve comes with a Weierstrass equation, and we find a 
quadratic twist of the base change of this curve to F which has conductor 9T. Since every 
elliptic curve over Q is modular, the base change and its quadratic twist are also modular, 
and since we exhaust the space this way we are done. 

For example, when 9T = (8) we have dimS l 2(8) new = 1 and find the form / with Hecke 
eigenvalues a^ w+2 ) = -2, a (3) = 2, a (w+3) = a^ w+7) = -4, a (lu _ 5) = a (w+4 ) = 4, and so on. In 
the tables of Cremona, we locate the form 

g(q) = q - 2q 3 - 2q 7 + q 9 - Aq 11 - 4g 13 - 4g 19 + ■■■ E S 2 (200) 

associated to an elliptic curve A g with label 200b and equation y 2 = x 3 + x 2 — 3x — 2. We 
find that the quadratic twist of (the base change of) A g by —(w + 2) = — Vo given by the 
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equation y 



x 3 + 



w 



l)x 2 — wx has conductor (8), so we conclude that / is the base 



change of g. Note that this example is not covered by the known cases of the Eichler-Shimura 
construction, since the level (8) = (2) 3 is squarefull. 

Example 6.5. Now consider the prime level 9? = (3w + 7) | 61, then S , 2(9 r t) new is a 2- 
dimensional irreducible Hecke module arising from the Jacobian of the Shimura curve Xq (1) 
where B is the quaternion algebra ramified at one of the infinite places and 9T. This follows 



from Examples |4.4| and 5.4 



One is naturally led to ask, can one identify the genus 2 Shimura curve X 
B of discriminant 9T? 

Consider the hyperelliptic curve C : y 2 + Q(x)y = P(x) over F with 

-wx 4 + (w- l)x 3 + (5w + A)x 2 + (6w + A)x + 2w + 1 



X B (1) with 



P{x) = - 

Q(x) = x 3 + (w — l)x 2 + wx + 1. 

This curve was obtained via specialization of the Brumer-Hashimoto family [U |3T] of curves 
(see also Wilson [73]) whose Jacobian has real multiplication by Q(a/5), so in particular so 
does the Jacobian J = Jac(C) of C. The discriminant of C is disc(C) = 9T 2 . One can show 
that J is modular using the theorem of Skinner and Wiles Theorem A] and the fact that 
J has a torsion point of order 31. Since C is hyperelliptic of genus 2, we compute that the 
conductor of (the abelian surface) J is 9T 2 and that its level is 91. Since J is modular of level 
91, it corresponds to the unique Hecke constituent of that level. It follows then by Faltings' 
isogeny theorem that J and the Jacobian of the Shimura curve X are isogenous. 



These examples exhaust all modular abelian varieties of GL 2 -type over 
01 of norm N{SR) < 100 (up to isogeny). 



with level 



We conclude with some (further) cases where Conjectures 6.1 and 6.2 are not known. 



Example 6.6. The first case of squarefull level with a form that is not a base change is level 
0T = (w + 36) = (w + 3) 2 so iV(9T) = 121. We have dim S2(w + 36) = 1 and the corresponding 
newform / has the following eigenvalues (where t 2 = 3): 

p (2) (w + 2) (3) (w + S) (w-A) (w-5) (w + 4) (w + 5) (to - 6) 
Np A 5 9 11 11 19 19 29 29 

a p (f) | t -t -1 -2t At -2 -5t -3 

In this case, the Eichler-Shimura conjecture predicts the existence of an abelian surface 
defined over F = Q(\/5) with good reduction away from p = (w + 3) and with real multipli- 
cation by Q(v3). The form / is not a base change from Q since a p (f) ^ ctp(f) in general. 
Can one find an explicit genus 2 curve over F, analogous to the previous example, with the 
L-function of its Jacobian given by the above Frobenius data? 

Example 6.7. We next consider a situation when we are successful in establishing the corre- 
spondence in a case which is not covered by known results. Consider 9T = (17w— 8) = (w+A) 2 
so Nfil) = 361. There exist forms f,g G S^^T) with Hecke eigenvalues as follows: 



p 


(2) 


(w + 2) 


(3) 


(w + 3) 


(w - 4) (w- 5) (w + 4) 


(w + 5) (w — 6) 


Np 


4 


5 


9 


11 


11 19 19 


29 29 


<h(f) 


2 


-3 


1 


3 


3-10 


3 -6 




-2 


-3 


-1 


-3 


-3 -1 


3 6 



2!) 



The form g is a quadratic twist of / by w(w + 4), and the forms f,g are not base changes 
from Q. 

We use the method of Cremona and Lingham [TU] — without an attempt to be exhaustive — 
to find an elliptic curve E over F with good reduction away from p = (w + 4). We find the 
curve 

E : y 2 + (w + l)y = x 3 + wx 2 — x — w 

(which could also be found by a naive search) with conductor 9T. We verify that #E(F p ) = 
Np + 1 — a p (f) for all primes p up to the limit of the computation. We prove that E is 
modular using the fact that E has a 3-isogeny, with kernel defined over F' = F(^w(w + 4)) 
generated by the points with x = w — 1, so we can apply the theorem of Skinner and Wiles. 
Therefore indeed E corresponds to /. 

We now turn to the existence of elliptic curves with everywhere good reduction over totally 
real fields. The Shimura-Taniyama conjecture predicts that such curves arise from cusp forms 
of level (1). 

Example 6.8. Let F = Q(A) = Q(Ci3) + be the totally real subfield of the cyclotomic field 
Q(Ci 3 ) with A = Cis + Cil 1 - 

The space 62(1) over F has dimension 1, represented by the cusp form /. Conjecture |6.1 
predicts the existence of an elliptic curve Af over F with everywhere good reduction. We 
note that Af must be a Q-curve since the twist f a of / by an element a G Gal(F/Q) is 
again a form of level 1 and so (since dimS^l) = 1) f a = f and consequently Af must be 
isogenous its Galois conjugates. 

Indeed, let 

A : y 2 + a\xy + a^y = x 3 + (I2X 2 + a^x + 

be the elliptic curve with the following coefficients: 

ai = X 2 + A + 1, 

a 2 = -A 3 + A 2 + A, 

a 3 = A 5 + A 3 + A, 

a 4 = A 5 - A 4 - 11A 3 - A 2 + 3A, 

a 6 = 13A 5 - 31A 4 - 48A 3 + 52A 2 + 33A - 10. 

We verify that A has everywhere good reduction and that #A(F) toIS = 19 with A(F) t0 rs 
generated by the point 

(A 5 - A 4 - 4A 3 + 3A 2 + 2A - 1, -A 5 + A 4 + A 2 ). 

Therefore, by Skinner and Wiles [61, Theorem 5] applied to the 3-adic representation at- 
tached to A or [62| Theorem A] to the 19-adic one, we conclude that A is modular, and thus 
A is our abelian variety Af. 

The elliptic curve A was obtained from a curve over Q(\/l3) C F as follows. There is a 
2-dimensional Hecke constituent g in the space of classical modular forms S l 2(169) new over Q 
given by 

g(q) =q + V^q 2 + 2q 3 + q 4 - V3q 5 + ...; 
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let A g be the associated abelian surface, defined over Q. Roberts and Washington [19] showed 
that A q twisted by is isomorphic over F to E x E a , where 



9 o 3483V13 - 12285 , r - 
E:y 2 = x 3 + ^— + (74844^13 - 270270), 



a Q-curve of conductor (13) and a is the nontrivial automorphism of Q(vL3)- 

The abelian variety A g is related to the elliptic curve A as follows. Roberts and Washington 
show that Ji(13) and A g are isogenous over Q. The abelian surface Ji(13) is isogenous over 
F to A x A, since Ji(13) obtains good reduction over F (see Mazur and Wiles [HH Chapter 
3.2, Proposition 2] and Schoof [52]). Indeed, Serre showed that the abelian surface Ji(13) 
has a rational point of order 19 and Mazur and Tate showed that Ji(13) twisted by Q(V^13) 
is a product of two elliptic curves. We thereby determine an equation for A by showing that 
it is a quadratic twist of E by some unit in F. We would like to thank Elkies and Watkins 
for pointing us to the rich history of this curve. 

The next example gives the list of all modular elliptic curves (up to isogeny) with every- 
where good reduction over a real quadratic field F of discriminant < 1000 and having strict 
class number one. 



Example 6.9. Let < D < 1000 be a fundamental discriminant such that F = Q(v_D) 
has strict class number one. We have 7Lp = Z[u>] with w = VT), (1 + \/D)/2, according 
as D = 0, 1 (mod 4). Let A be a modular elliptic curve over F with everywhere good 
reduction. Then A corresponds to a cusp form / of level 1 with rational Fourier coefficients. 
Computing the spaces 62(1), we find that either A is a Q-curve (so A is isogenous to its 
Galois conjugate), or D G {509,853,929,997} and the isogeny class of A is represented by 
the curve with coefficients as below, up to Galois conjugacy: 



D 


ai 


a 2 




a i 


a 6 


509 


-1 


2 + 2w 


— w 


162 + 3w 


71 + 34w 


853 


w 


— 1 — w 





-5921 + 429w 


145066 - 9437w 


929 


w 


1 — w 


1 + 20w 


-1738 - 82w 


-11808 -201w 


997 





w 


1 


-289 - 2Aw 


-2334 - lUw 



The Q-curves over F are listed in Cremona [8]: for example, the curve y 2 + xy + uy = x 3 
with u = (5 + v / 29)/2 is the smallest such and was discovered by Tate [54] . 

The curve for discriminant D = 509 was discovered by Pinch [52, p. 415]: it was the first 
example of an elliptic curve having everywhere good reduction which is not a Q-curve, and 
it was proven to be modular by Socrates and Whitehouse [60] using the method of Faltings 
and Serre. The curve for D = 853 was discovered by Cremona and Watkins independently; 
the one for D = 929 by Elkies; and the one for D = 997 by Cremona. Their modularity 
can be established by applying the theorem of Skinner and Wiles [61, Theorem 5] to the 
associated 3-adic representation. 

As a third and final application, we consider the existence of number fields with small 
ramification — these are linked to Hilbert modular forms via their associated Galois repre- 
sentations. 

In the late 1990s, Gross proposed the following conjecture. 
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Conjecture 6.10. For every prime p, there exists a nonsolvable Galois number field K 
ramified only at p. 



Conjecture 6.10 is true for p > 11 by an argument due to Serre and Swinnerton-Dyer 
For k = 12, 16, 18,20,22 or 26, let A k e S/-(l) be the unique newform of level 1 and 
weight k. For each such k and every prime p, there is an associated Galois representation 
p kp : Gal(Q/Q) — > GL 2 (F P ) ramified only at p. In particular, the field K fixed by the kernel 
of pk tP is a number field ramified only at p with Galois group Gal(_fT/Q) = imgp fcp . 

We say that the prime p is exceptional for k if p k is not surjective. The "open image" 
theorem of Serre [55] implies that the set of primes where p k is exceptional is finite. Serre 
and Swinnerton-Dyer produced the following table of exceptional primes: 



_k_ 
12 
16 
18 
20 
22 
26 



Exceptional p 
2,3,5,7,23,691 
2,3,5,7,11,31,59,3617 
2,3,5,7,11,13,43867 
2,3,5,7,11,13,283,617 
2,3,5,7,13,17, 131,593 
2,3,5,7,11,17,19,657931 



(The large primes that occur are divisors of denominators of Bernoulli numbers.) It follows 



that for every p > 11 that Conjecture 6.10 is true. 

Conversely, it is a consequence of the proof of Serre's conjecture by Khare and Winten- 
berger [39] together with standard level lowering arguments that if p < 7 then any odd 
representation of the absolute Galois group Gal(Q/Q) to GL 2 (F P ) which is ramified only at 
p is necessarily reducible (and thus solvable). 

To address Gross' conjecture, we instead look at the residual representations associated 
to Hilbert cusp forms of parallel weight 2 and level 1 over fields F ramified only at a prime 
p < 7. This idea was first pursued by the first author [16] following a suggestion of Gross, 
who found such a form with F = Q(C32) + , the totally real subfield of the cyclotomic field 
Q(^32), settling the case p = 2. This line of inquiry was followed further by the authors with 
Matthew Greenberg [T7] . settling the cases p = 3,5. (The case p = 7 was recently settled by 
Dieulefait [2U] by considering the mod 7 Galois representation attached to a genus 2 Siegel 
cusp form of level 1 and weight 28. Therefore Gross' conjecture is now a theorem.) 

We sketch below the resolution of the case p = 5. We take the base field F = Q(6) to be 
the subfield of Q(C2s) of degree 5, where 6 5 + 56 4 — 256 2 — 256 — 5 = 0. The field F has strict 
class number 1. Let E be the elliptic curve over F with j-invariant given by 

(5 • 7 n )j = - 168635243727774766 4 - 885403699379835886 3 + 112479 1466055321 5b 2 
+ 4643993605154835726 + 353505866738383680 

and minimal conductor 91 Then 9? = p 5 p 7 , where p 5 = ((-26 3 - 126 2 + 316 + 25)/7) is the 
unique prime above 5 and pj = ((— 26 4 — 96 3 + 86 2 + 536 + 6)/7) is one of the five primes 
above 7. Roberts [50] showed that the mod 5 Galois representation 

p E ^ : Gal(F/ F) -> End(£[5]) S GL 2 (F P5 ) = GL 2 (F 5 ) 

has projective image PGL^Fs) and is ramified only at the prime p5 (and not P7). (We will 
verify below that E has indeed the right conductor.) The representation p E5 and its Galois 
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conjugates gives an extension K of Q with Galois group 

Gal(AT/Q) = PSL 2 (F 5 ) 5 x Z/10Z. 
By work of Skinner and Wiles (as in Example 6.3[ ), we prove that the 3-adic representation 



Pe,3 associated to E is modular, and hence E itself is modular. Since E is modular and F 
has odd degree, E is uniformized by a Shimura curve. Namely, let B be the quaternion 
algebra over F ramified at four of the five real places. Let X§ (P5P7) be the Shimura curve 
associated to an Eichler order of level p 5 p 7 contained in a maximal order O, and let Jq (P5P7) 
be the Jacobian of X^(p 5 p 7 ). We compute that Jq (psp7) new has dimension 203. We then 
find the unique Hilbert newform fp of parallel weight 2 and level P5P7 with integer Fourier 
coefficients which corresponds to E. 

The elliptic curve E found by Roberts [50] was obtained from our computations of Hilbert 
modular forms at level ps over F. The space S2(p5) new has 2 Hecke constituents of dimension 
10 and 20, respectively. Let / be a newform in the 20-dimensional constituent. Let Tf be 
the restriction of the Hecke algebra T = Z[T P ] P to this constituent: this is the constituent 
which yields the Galois representation obtained in the 5-torsion of Roberts' curve. Let 
E f = T f ® Q = Q(a p (f)) be the field of Fourier coefficients of / and let A = Aak(E f ). By 
direct calculations, we see that Kf = E^ is the (totally real) quartic field generated by a 
root of the polynomial x 4 + 2x 3 — 75x 2 — 112x + 816. 

The Galois group Gal(F/Q) acts on T via its action on ideals of Tp, namely &(T p ) = T CT ( p ). 
This action preserves the decomposition of T into Hecke-irreducible components. By work 
of Shimura, the action of A on Fourier expansions preserves the Hecke constituents. In 
particular, both these actions preserve Tf and hence Kf and must be compatible. Therefore, 
for each a G Gal(F/Q), there is a unique r = r(er) £ A such that, for all prime ideals p C Tp, 
we have 

<Mp)(/) = r(a p (f)). 
The map r thus yields a homomorphism 

Gal(F/Q)^A 
(611) a^r(a) 

By direct calculation, we show that this map is an isomorphism. Since Gal(F/Q) is abelian, 
the field Ef must be a ray class field over Kf] in fact, we compute that it has conductor equal 
to a prime above 71 which splits completely in Ef. The ideal 5Z#. factors as 5T^ f = *p 2( P'. 
The prime *P splits completely in Ef, and the primes above it in Ef are permuted by 
Gdi(Ef/K f ) = Aut(E f ); the prime *p ; is inert. 
We have a Galois representation 

p f : Gal(F/F) GL 2 (T / ®F 5 ) 

with ^(^(Frobp)) = a p (f) (mod 5) and det(p / (Frob p )) = A^p (mod 5) for all p \ 5. Let 

vcif 1 , i — 1, . . . , 5 be the five maximal ideals above (the image of) They give rise to the 5 
residual representations 



pf : Gal(F/F) -> GL 2 (F 



5 j 



Let Li be the fixed field of ker(pf), and L the compositum of the L 4 . Then L is a Galois 
extension of F and is ramified only at 5. Since F is Galois and ramified only at 5, and by the 
above Gal(F/Q) permutes the fields Lj so preserves L, we see that L is a Galois extension 
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of Q ramified only at 5. By a result of Shepherd-Barron and Taylor [56], each extension Li 
can be realized in the 5-torsion of an elliptic curve Ei/F. 

Recall the projective representation Pp E 5 from Roberts' elliptic curve E is surjective, 
ramified at p 5 but unramified at P7. Therefore, the level p 5 p 7 is a nonoptimal level for p E5 ; 

thus, by Mazur's Principle [37], we have ~Pe<o — Pf f° r some 1 < i < 5. In other words, the 
extension constructed by Roberts is isomorphic to our field L. 

Roberts |50j has given an explicit equation for the number field L (obtained from the 
5-division polynomial of the elliptic curve E): the field L is the the splitting field of the 
polynomial 

a; 25 - 25a; 22 + 25a; 21 + 110a; 20 - 625a; 19 + 1250a; 18 - 3625a; 17 
+ 21750a; 16 - 57200a; 15 + 112500a; 14 - 240625a; 13 + 448125a; 12 

- 1126250a; 11 + 1744825a; 10 - 1006875a; 9 - 705000a; 8 + 4269125a; 7 

- 3551000a; 6 + 949625a; 5 - 792500a; 4 + 1303750a; 3 - 899750a; 2 
+ 291625a; - 36535. 



To conclude, we consider a question that touches on each of the above three subjects. We 
reconsider Gross' conjecture (Conjecture 6.10) in the case p = 2 [16J. The nonsolvable field 
which is ramified only at 2 arises from the Galois representation associated to a constituent 
eigenform / £ £2(1) in a 16-dimensional subspace of the space of Hilbert cusp forms of 
parallel weight 2 and level 1 over F = Q(C32) + - Let Ef = Q(a p ) be the field of Fourier 
coefficients of / and let A = Aut(Ef). Let Kf be the fixed field of A so that Gal(Ef / Kf) = 
A. The map (6.11) in this context is again an isomorphism, and so Ef is abelian over Kf 



and 8 = [F : Q] divides [Ef : Q] = 16. By direct calculations, we show that Kf = 

The Eichler-Shimura construction (Conjecture |6. 2 ) predicts the existence a 16-dimensional 
abelian variety Af defined over F with everywhere good reduction and real multiplication 
by Ef associated to /. More should be true, as communicated to us by Gross (private 
communication). In fact, Af should descend to an abelian variety of dimension 16 over Q, 
and we should have L(Af,s) = L(f,s)L(f a ,s), where a £ Gal(-E//Q) is any element that 
restricts to the nontrivial element of Kf = Q(a/5). The endomorphisms of Af over Q should 
be the ring of integers of Kf, and over Q, the variety Af would have bad reduction only at 
the prime 2; the nonsolvable extension would then arise as its 2-division field. The conductor 
N = 2 124 = c? 4 of Af over Q can be computed from the functional equation of L(f, s), where 
we note d = 2 31 is the discriminant of F. We note that Af is not of GL/2-type over Q itself 
as it would be modular by the proof of Serre's conjecture. 

Although one knows that the associated Galois representation exists by work of Taylor 
[61)] . as Gross says, 'With such nice properties, it's a shame that we can't even prove that 
the abelian variety Af exists! That's an advantage you have when F has odd degree." 



7. Adelic quaternionic modular forms 

In this section, we begin again, and we revisit the definition of Hilbert and quaternionic 
modular forms allowing F to have arbitrary class number; we refer to Hida [31] as a reference 
for this section. 
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We renew our notation. Let F be a totally real field of degree n = [F : Q] with ring of 
integers "Lp. Let B be a quaternion algebra over F of discriminant D. Let v\, . . . , v n be the 
real places of F (abbreviating as before Xj = Vi(x) for x G F), and suppose that B is split 
at Vi, . . . , v r and ramified at iv+i, ■ ■ ■ ,v n , i.e. 

(7.1) BHB 0O = B® (J l4M 2 (l) r xE n - r . 

Let ^ denote the ith projection B — > M 2 (R) and i = (ti, . . . , t r ). Let 

F+ = {s G F : 2j > for all i} 

be the group of totally positive elements of F and let Z£ + = Z F fl F_*. 

Let Oo(l) Q B be a maximal order. Let OT be an ideal of T,p coprime to D and let 
O = e> (0?) C O (l) be an Eichler order of level 9t 

With a view towards generalizations, rather than viewing modular forms as functions on 
(a Cartesian power of) the upper half-plane which transform in a certain way, we now view 
them instead more canonically as functions on B^. Let = C \ R be the union of the 
upper and lower half-planes. Via the embeddings Vi,...,v r , corresponding to the first r 



factors in (7.1), the group B^ acts on ("H ± ) r on the right transitively with the stabilizer of 
(v^T, . . . ,7^1) e U r being 

= (M x S0 2 (M)) r x (H x ) n ~ r . 

Therefore we can identify 
(7.2) 



B^/K^ -> (^Y 



From this perspective, it is natural to consider the other (nonarchimedean) places of F at 
the same time. Let 

/ 

Z = hmZ/nZ = Yl % P 

n p 

(where ' denotes the restricted direct product) and let ~ denote tensor with Z over Z. We 
will define modular forms on B as analytic functions on B^ x B x which are invariant on 
the left by B x and transform by x O x in a specified way. 

We must first define the codomain of these functions to obtain forms of arbitrary weight. 
Let k = (ki, . . . , k n ) G (Z>2) n and suppose that the ki have the same parity; we call k a 
weight. Let 

ko = maxki, rrii = (&o — ki)/2, and W{ = ki — 2. 

i 

For an integer w > 0, let P w = P W (C) be the subspace of C[x, y] consisting of homogeneous 
polynomials of degree w. For 7 G GL 2 (C), let 7 be the adjoint of 7, so that 77 = 77 = det 7. 
Define a right action of GL 2 (C) on P W (C) by 

(q ■ 7) (2, y) = q((x y)j) = q(dx - cy, -bx + ay) 

for 7 = ^ G GL 2 (C) and q G P W (C). For m G Z, GL 2 (C) also acts on F„,(C) via 

the character 7 1— >■ (det7) m . By twisting the above action by this character, we get a right 
GL 2 (C)-module denoted by P w (m)(C). Define the right GL 2 (C) n - r -module 

(7.3) W k (C) = W(C) = P Wr+1 (m r+1 )(C) ® • • • ® F U)n (m n )(C). 
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(By convention, if r = n then we set Wfc(C) = C.) For the ramified real places tv+i, . . . ,v n 
of F, we choose splittings 

ii:B4B® f C^M 2 (C). 

We abbreviate as above 7$ = ^(7) for 7 G B. Then Wk(C) becomes a right .B x -module via 
7 1 — ^ (jr+i: ■ ■ ■ )7n) £ GL 2 (C) n_r . From now on, Wjt(C) will be endowed with this action, 
which we denote by x ^ i 7 for x G W&(C) and 7 65 x . One may identify W&(C) with the 
subspace of the algebra C[x r+ i, y r +i, ■ ■ ■ , ^n, J/n] consisting of those polynomials q which are 
homogeneous in (xi,yi) of degree Wi but with a twisted action. 

We consider the space of functions : B^ x B x — > Wk(C), with a right action of x B x 
defined by 



(7.4) (*i*(ic,3))Gfc5) = fn ( ^;ff-i j 

where recall j(/tj, a/— T) = Ci\f— 1 + <ij G C if «j = ^ l . (The presence of the inverses is 

forced as we want a right action on functions via multiplication of the argument on the right. 
This almost extends to a right action of x B x , except that j(gh, z) = j(g, hz)j(h,z) 7^ 
j(g,z)j(h,z) unless h fixes z.) 

Definition 7.5. A (quaternionic) modular form of weight k and level 91 for B is an analytic 
function 

: 5 X x 5 X ^ W k (C) 
such that for all (g,a) G B^ x B x we have: 

(i) (0 \k (k, u))(g, a) = 0(g, a) for all k G and u G O x ; and 

(ii) 0(75-, 7S) = 0(^(, a) /or a// 7 G 5 X . 



In other words, Definition |7.5| says that a quaternionic modular form of weight k and level 
is an analytic function which is 5 x -invariant on the left and (K^ x (9 x )-equivariant on 
the right under the action (7.4). In particular, we have 

(7.6) <f>(9,aru) = (<p\ k (l,u))(g,au) = (j)(g,a) 
for all u G O x and 

(7.7) fan, a) = (0 \ h (k, 1))( 9 k, a) = (f[ ! x 1 fa, 

We denote by Mf (01) the space of quaternionic modular forms of weight k and level 9T 
for B. 

Let be a quaternionic modular form. For (z,aO x ) G ('H ± ) r x B x /O x , choose g G B x 
such that g{\/— 1, • • • , \/— 1) = 2 and define 

/ : (H ± ) r x £?7<5 X -4 W fc (C) 
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OO 



The map / in (7.8) is well-defined by (7.6) and (7.7): we have 
(det gi Ki) m ^- 1 



n 



\ j(9i K i 



-l)ki 



<f>(gK,auY 9Kr 



n 



(det g { 



mi+ki — l 



<t>{g, S) 



using the fact that j{giK h a/-1) = j(g h y/^l)j{K h v 7 -!)- 

The identity 7.5(H) is translated as follows. Let 7 G B x . Then g(y/—l, 

— 1) = 72;, so 



2; IS 



equivalent to (jg) (y— T 
/( 7 2:,7a6 >< 

(7.9) 



n 



■7'(7i^» V-l) 
J- (det 7^, 



(79)" 



n 



1 (det 7^ 



mi+ki-l 



n 



1 (det fl r i )'ni+*<-l 



( /)(^,a) i 



n 



j(t»>«. 



\ ( det 7i) 



where now Zi) = CiZ { + dj G C if j { 



Ci (ij 



for all 7 G i? x and we have the relation 



j(^yS, z) = j(7, 8z)j(5, z) for all z G "H and 0, /1 G 5 X . Accordingly, we define a right action 
of B x on the space of functions in (7.8) by 



(7.10) 
Then 



(/U 7 )(z,30> 



n 

v i=i 



(det 7i 



imi+fej— 1 



j(7i>^ 



/( 7 z,7«0 X ) 7 - 



(/| fc 7)(^«0 X ) = /(^«0 X ). 

Note that the central subgroup F x C 5 X acts by (/ | fc a)(z, q?C x ) = N F / Q (a) k °~ 2 f(z, aaO x ) 
for a e F x . 

The C-vector space of modular forms of weight k and level 91 on B is finite-dimensional 
and is denoted Mf (91). 

Lemma 7.11. There is a bijection between M^(D1) and £/ie space of functions f : ("H ± ) r x 
B x jO x Wfc(C) £/iai are holomorphic in the first variable and locally constant in the second 
one and such that 

f\kl = f for all <ye B x . 



From now on, we will only work with modular forms / as presented in Lemma 7.11 



We define the quaternionic Shimura variety of level 01 associated to B as the double coset 

X B (D1)(C) = B X \{B X JK^ x B X /6 X ) = B x \((^ ± ) r x 5 x /a x ); 

the set Xq(%1)(C) can be equipped the structure of a complex (possibly disconnected) Rie- 
mannian manifold of dimension r. 

Example 7.12. We recover first the definition of classical modular forms when F = Q and 
B = M2(Q). For simplicity, we take N = 1. In this case, r = 1 so Wk(C) = C and m — 0. 
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The action (7.10) is simply 



(/ \ h 7 )(*, «a x ) = (d f t) " x f{ iz , 7 «a x ). 



We take the definition (|7.8|) as our starting point. The element ( ^ ^ ] G B x = GL- 



identifies the upper and lower half-planes, so a modular form / : 'H ± x GL 2 (Q)/ GL 2 (Z) — > C 
is determined by its restriction / : % x GL 2 (Q)/T (A) — >■ C and the subgroup of GL 2 (Q) 
which preserves % is exactly GL^(Q). 

We wish to recover the classical action using this new action of GL 2 (Q), so we are led 
to consider the double coset GL^ (Q)\GL 2 (Q)/ GL 2 (Z). An element of this double coset is 
specified by an element a G GL 2 (Q) up to right-multiplication by GL 2 (Z), i.e. a Z- lattice 
in Q 2 , i.e. A G Lat(Q 2 ), specified by the rows of a. But the map Lat(Q 2 ) ->■ Lat(Q 2 ) by 
A i— > A is a bijection, with inverse A^AnQ 2 . And since GL^ (Q) acts transitively on the 
left on the set of lattices Lat(Q 2 ), we conclude that # GL^ (Q)\GL 2 (Q)/ GL 2 (Z) = 1. 

It follows that / is uniquely specified by the function f(z, O x ) for z G %, which by abuse 
of notation we write simply /:%—»• C. The stabilizer of GL 2 (Q) acting on GL 2 (Z) by 
multiplication on the left is GL 2 (Z), so we recover the condition 

3K7, z ) 

for all 7 G GL 2 (Z), which is exactly the definition given in Section 1. 

The interested reader can modify this argument for N > 1; alternatively, we give a general 
derivation below. 

We now define cusp forms. If B = M 2 (F), then we are in the situation of Hilbert modular 
forms (but over a field with arbitrary class number): so we define a cusp form to be a form 
/ such that f(z) — > whenever z tends to a cusp P 1 (F) <^-> P 1 (IR) ?1 . Otherwise, the Shimura 
variety X^(0 r t)(C) is compact. If B is indefinite, so < r < n, then there are no cusps, and 
we define the space of cusp forms to be Sj?(yi) = M^(9T). Finally, suppose B is definite; 
then r = 0. If k = (2, . . . , 2), then we define E% {% to be the space of those / G Mjf {% 
such that / factors through nrd : B x — > F x ; otherwise, we set E^i^l) = 0. Then there is an 
orthogonal decomposition Afjf = Sj* © E%(% and we call Sj* (9?) the space of cusp 
forms for Z?. 

The spaces M^(OT) and S^iVl) come equipped with an action of pairwise commuting 
diagonalizable Hecke operators T n indexed by the nonzero ideals n of Z^, defined as follows. 
Given / G Sj?(yi) and n E B x , we define a Hecke operator associated to 7? as follows: we 
write 



(7.13) 6 x 7r6 x = I I (3 



x ~ 
7T,- 



and let 

(7.i4) (T^)(^,sa x ) = ^/(z.aTfr^: 
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(Again, although it may seem unnatural, the choice of inverse here is to make the definitions 
agree with the classical case.) 

For a prime p of with p \ 2)91, we denote by T p the Hecke operator 7~ where tt G B x 

is such that tt v = 1 for v 7^ p and n p = ^ jj G <E>z F Z^p = M 2 (Z F) p) where p 6 Z FiP is 

a uniformizer at p. 

Equivalently, for a prime p and p G Zj? such that pZ ^ fl Z^ = p , we define 

9(p) = O x \{tt G O : nrd(Tr) = p}. 
a set of cardinality iVp + 1, and define 

(7.15) (T p f)(z,aO x ) = fiz^O*) 

?ee(p) 

where we have implicitly chosen representatives tt G O for the orbits in 0(p). For an ideal 
n of Tip, the operator T n is defined analogously. 

We say that a cusp form / is a newform if it is an eigenvector of the Hecke operators 
which does not belong to M k (Tl) for Wl | 9t. 

To unpack this definition further, and to relate this definition with the definitions given 
previously, we investigate the structure of the Shimura variety 

X B (9t)(C) = B X \(B^ /K 00 x B X /6 X ) = B x \({U ± ) r x B x /d x ). 

By Eichler's theorem of norms [66, Theoreme III.4.1], we have 

nrd( J B x ) = F x +) = {a G F x : u f (o) > for z = r + 1, . . . , n}, 

i.e. the norms from B x consists of the subgroup of elements of F which are positive at all 
real places which are ramified in B. In particular, B x /B x = (Z/2Z) r , where 

B x = { 7 G B : nrd( 7 ) G F x } 

is the subgroup of B x whose elements have totally positive reduced norm. 
The group B+ acts on ~H r , therefore we may identify 

X B (9T)(C) = B x \(H r x B X /3 X ) 

and a modular form on ('H ± ) r x B x jO x can be uniquely recovered from its restriction to 
W x B x jO x . Now we have a natural (continuous) projection map 

X B (9T)(C) ^B X \B X /3 X . 
The reduced norm gives a surjective map 

(7.16) nrd : B X \B X jO x F X \F X /Z x ^ Cl + Z F . 

where Cl + Z F denotes the strict class group of Z F , i.e. the ray class group of 7L F with modulus 
equal to the product of all real (infinite) places of F. The theorem of strong approximation 



[661 Theoreme III. 4. 3] implies that (7.16) is a bijection if B is indefinite. So our description 
will accordingly depend on whether B is indefinite or definite. 

First, suppose that B is indefinite. Then space X B (9T)(C) is the disjoint union of connected 
Riemannian manifolds indexed by Cl + Zi?, which we identify explicitly as follows. Let the 
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ideals aCZ f form a set of representatives for Cl + Z F , and let a G Z F be such that a Z F fl 
Zjr = o. (For the trivial class a = Zp, we choose a = 1). By (7.16), there exists a G -B x such 
that nrd(a) = a. We let O a = ada^DB so that C (1) = £>, and we put T a = O x + = 3 x nB x . 
Then we have 



(7.17) 



x B (*n)(C)= □ Bi(H r xad*)^ □ r \?r, 

[a]sCi+(z F ) [a]eci+(z F ) 
where the last identification is obtained via the bijection 

(7.18) + XV ^ ' X 

B x (z,aO x ) m- z 

Now let / G Afjf (91), so that / : (H ± ) r x £ x /^ x ~> W k (C) satisfies / | fe7 = / for all 
7 G B x . Let M^(0T, a) be the space of functions f a \W^- W k (C) such that /„ \ k 7 = /„ for 
all 7 G r„, where we define 



n 



(det 7i) 



mi+fcj— 1 



for 7 G -B x . Then by (7.18), the map 
(7.19) 
where 



[o]eci+(z F ) 
/ H- (/ a ) 



/a : -> W fc (C) 
/„(*) = /(z,a0> 



is an isomorphism. 



We now explain how the Hecke module structure on the left-hand side of (7.19), defined in 



(7.14)-(7.15), is carried over to the right-hand side. We follow Shimura [57, Section 2]. We 
consider the action on the summand corresponding to [a] G Q + (Zp). Extending the notation 
above, among the representatives chosen, let b be such that [b] = [dp -1 ], let bZ F fl 7L F = b, 
and let /3 G B x be such that nrd(/3) = b. 
By definition, 

(T p f) a (z) = (T p f)(z,ad x ) = f{z,™- X 6*)- 

TfeeO) 

Let 7? G @(p). Then by strong approximation, we have 

(7.20) mr^o/r 1 n B = w- 1 O b 

with w G B x , since this lattice has reduced norm [ap _1 (ap _1 ) _1 ] = [(1)]. Therefore, there 
exists u G O such that a7r _1 M/3 _1 = w^ 1 whence 

(7.21) (T p f) a (z) f{z^ l 6 x ) = Y,f( z ^- l PO x Y 

7r60(p) 
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the second sum runs over a choice of w as in equation (7.20) corresponding to each tt G 0(p). 
This latter sum can be identified with a sum over values of A as follows. We have 

f(z, w-^O*) = (/ \ k w)(z, w- x W) 

i= l Jv^ii z i) % J 

The first equality follows from the 5 x -invariance of / and the others by definition of the 
slash operators. Putting these together, we have 

(7.22) (T p f) a (z) = J2(h\k^)- 

(The naturality of this definition explains the choice of inverses above.) 

This adelic calculation can be made global as follows. Let I a = aO fl B and h = /30 fl B. 
For n G 0(p), we have 

mT^/r 1 HB = w~ x O b 

hence 

(7.23) O b zu = pdna- 1 = (pd)a~ l (aOcT^wrcT 1 nB = hi; 1 J. 

The elements w thus obtained are characterized by their norms (in the right lattice), as with 
the Hecke operators defined previously (5.1): we analogously define 

e(p) a ,b = r b \{zue hi; 1 n B$ : nrd^/^p = (nrd(w))} 

= T b \ {w G hi; 1 n B* : nrd(tt;)b = ap}. 

Then for A e Mjf (*Tl, b), we have T p A e M%{% a) and 



6 |fc ^ 

where [b] = [ap -1 ]. 

Example 7.25. If F has strict class number 1, then I a = I b = O = O a so O b w a = OiT a as in 
Section 5. 



We note that the isomorphism (7.19) preserves the subspace of cusp forms in a way that 
is compatible with the Hecke action, so we have a decomposition 

[o]eci+(z F ) 



Example 7.26. Let B = M 2 (F), and let O = O {% C O (l) = M 2 (Z F ). Then we may take 
a 0" 
1 



a — ( ? ? 1 G GL 2 (-F), and so we find simply that 



Let 



a = sm^z^s- 1 n B = (1^ *^ = O {% a) 
r (%a) = O (m,a) x + = | 7 = r ^ G o (OT,o) : det 7 G Z* + J . 



41 



Then 

X ( f(OT)(C)= □ r (%a)\H n 

[a]eCl+(Z F ) 

is a disjoint union. 

A Hilbert modular form of weight k and level 91 is a tuple (/ a ) of holomorphic functions 
f a : T-L n — > C, indexed by Cl + Z,p, such that for all a we have 

(/a|*7)(*) =/«.(*) for all 7 Gr (0T,a) 

(with the extra assumption that / is holomorphic at the cusps if F = Q). Or, put another 
way, let Mfc(9T, a) be the set of holomorphic functions "H n — > C such that (/ \kj)(z) = f(z) 
for all 7 G r (9T, a); then 

M k (Vl) = 0M fc (9t, a). 

[■»] 

In particular, we recover the definitions in Section 2 when F has strict class number 1. 
A modular form / G Mfc(9T, a) admits a Fourier expansion 



f(z) = a + ^2 a M e " 



D 2iriTr(nz) 
Me(aa- 1 )+ 

analogous to ( |2.5[ ). We say that / G M fc (9T, a) is a cusp form if /(z) — > as z tends to any 
cusp. Letting Sjffll, a) be the space of such cusp forms, we have 

Sk(Vl) = @S k (% a). 

M 

Let / = (/„) G Sk(yi) be a Hilbert cusp form and let n C Z F be an ideal. Suppose that 
[n] = [aD -1 ] amongst the representatives chosen for Cl + (Z j p), and let v G be such that 
n = uad^ 1 . We define a n = i/ m a„(/ a ); the transformation rule implies that a n only depends 
on n and we call a n the Fourier coefficient of / at n. 

Now suppose that B is definite. Then the Shimura variety is simply 

X B (9t)(C) = B x \B x /O x — CIO 
and so is canonically identified with the set of right ideal classes of O. Note that the reduced 



norm map (7.16 ) here is the map nrd : CIO — > CI Z^ which is surjective but not a bijection, 
in general. A modular form / G M^ftl) is then just a map / : B x /O x — > Wfc(C) such that 
f\kl = / for all 7 G B x . Such a function is completely determined by its values on a set 
of representatives of CIO; moreover, given any right ideal / = aO fl B, the stabilizer of 
B x acting on aO by left multiplication is Ol{I) x = aO x a~ 1 fl B x . Therefore, there is an 
isomorphism of complex vector spaces given by 

M^{%^ ^ fc (C) r C) 

v ' ; I=a6nB 

f * (/(«)), 

where T{I) = SC^S" 1 n B x = O l {I) x and W k {C) T{I) is the r(J)-invariant subspace of 
W k (C). 
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Having now discussed both the definite and indefinite cases in turn, we return to a general 
quaternion algebra B. Let / G S 1 ^ (9T) ncw be a newform. A theorem of Shimura states that 
the coefficients a n are algebraic integers and Ef = Q({a n }) is a number field. The Hecke 
eigenvalues a n determine the L-series 

Hf,s)= = n f 1 - j^p-s + Np 2s + i- k0 ) nf 1 -^) 

ncz F pf*n v r r 7 p|m v r 7 

associated to / (defined for Res > 1). Moreover, associated to / is a Galois representation: 
for [ a prime of and the completion of Ef at [, there is an absolutely irreducible, 
totally odd Galois representation 

p u : Gal(F/F) GL 2 (£ /jt ) 

such that, for any prime p f [9t, we have 

Tr(p /it (Frob p )) = a p {f) and det(p / , [ (Frob p )) = A^ " 1 . 

The existence of this representation is due to work of Blasius-Rogawski [2], Carayol [6], 
Deligne [Lj, Saito [EE], Taylor [65], and Wiles 



The statement of the Jacquet-Langlands correspondence (3.9) reads the same in this more 
general context. 

Theorem 7.28 (Jacquet-Langlands). There is an injective map of Hecke modules 

whose image consists of those forms which are new at all primes dividing D. 
We are now ready to state the main general result of this article, generalizing the result 



of Theorem 2.8| to arbitrary class number and arbitrary weight. 



Theorem 7.29 (Dembele-Donnelly [IE], Voight [69]). There exists an algorithm which, given 
a totally real field F , a nonzero ideal 91 C 7L Fl and a weight k £ (Z> 2 )f F: ^, computes the 
space Sk{yi) of Hilbert cusp forms of level 91 over F as a Hecke module. 

The proof of this theorem is discussed in the next two sections. It falls again naturally 



into two methods, definite and indefinite, which overlap just as in Remark 3.10 

8. Definite method, arbitrary class number 

In this section, we return to the totally definite case but allow arbitrary class number. As 
explained above, the space Mjf (91) of modular forms of level 91 and weight k on B is the 
space of functions / : B x /O x ->• W k (C) such that / \ k j — f for all 7 G B x . 



We can use the identification (7.27) to compute the space Sjf (9t) as in the direct approach 
of Section 4, with the appropriate modifications. Let Ii, . . . , I H be a set of representatives 
for CI O such that nrd(Jj) is coprime to £>9t for all i. Let dti G O be such that a.iO D O = Ii, 



and let Tj = Oiih)* ■ Then dualizing the isomorphism (7.27), we have 

H 

M fc B (9r)^0H4(C) r \ 

i=l 

The Hecke module structure on this space is defined similarly as in Section 4, as the following 
example illustrates. 
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Example 8.1. Consider the totally real quartic field F = Q(w) where w 4 — 5w 2 — 2w + 1 = 0. 
Then F has discriminant 5744 = 2 4 359 and Galois group S4. We have Cl + "Lp = 2 (but 
C1Z F = 1). 

The quaternion algebra B = ^ — ^ — j is unramified at all finite places (and ramified at 

all real places). We compute a maximal order O and find that #C1(9 = 4. We compute 
the action of the Hecke operators as in (3.8): we identify the isomorphism classes of the 



iVp + 1 right ideals of norm p inside each right ideal J in a set of representatives for CI O. 
We compute, for example, that 

/0 




T, 



(w 3 — 4w— 1) 



2 

V 3 






2 
3 



4 







where N(w 3 — Aw — 1) = 4; note this matrix has a block form, corresponding to the fact that 



iw 



Aw 



1) represents the nontrivial class in CI Zp. Correspondingly, 

/6 2 \ 
5 12 



T, 



{w 2 —w—i) 












6 



with N(w 2 



w 



A) = 13 is a block scalar matrix, as 



4 

10/ 

2 



4) is trivial in C\ + Z*f- In 



this case, the space E 2 (l) of functions that factor through the reduced norm has dimension 
dim£ , 2 (l) = 2, so dimS , 2 (l) = 2, and we find that this space is irreducible as a Hecke module 
and so has a unique constituent /. 

We obtain the following table of Hecke eigenvalues: 



p 


(w 3 - Aw - I) 


(w-1) 


(w 2 -w- 2) 


(w 2 - 3) (w 2 - w - 4) 


(w 2 - 2) 


Np 


A 


5 


7 


13 13 


17 


<h(f) 





t 


-2t 


-t A 


3t 



Here t satisfies the polynomial t 2 — 6 = 0. As in Section 6, we predict the existence of an 
abelian variety over F with real multiplication by Q(v / 6) and everywhere good reduction. 



As in Section 4, the disadvantage of the approach used in Example 8.1 is that for each 
level 91, one must compute the set of ideal classes CIO = Cl(9o(91)- By working with a 
more complicated coefficient module we can work with ideal classes only with the maximal 
order Oq(1), as follows. 

Changing notation, now let h, . . . ,Ih be representatives for CI Oq(1), with h = # CI Oq(1), 
and let Jj = ctjC^l) H B. By strong approximation, we may assume that each nrd(Jj) is 
coprime to 2)91: indeed, we may assume each nrd(ij) is supported in any set S of primes 
that generate C1 + Z F . Let O (l)i = L (I t ) = a^o^Sr 1 n B be the left order of I t . Then 
O (l)i ® Zf Z F ,m = Oo(l) ®-l f %F,m- 

Let a for a G P 1 (Zi?/01) represent the Oo(l), O-connecting ideals of norm 01: that is, 
(3 a e 6 (1) and if J a = O {l)X n 5 then O fi (J a ) = O. Then the set {/;,/„};,„, where 
JjJ a = ai/Sa^ H -B, covers all isomorphism classes of right (9-ideals, but not necessarily 
uniquely: two such ideals IiJ a and IjJb are isomorphic if and only if i = j and there exists 
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7 G Oo(l)f such that jJ a = J&, comparing the elements ai(3 a ,aj[3b G 0o(l). The action of 
O (l) x can be equivalently given on the set of indices a G P 1 (Zjr/O r t): via the (reduction of 
a) splitting map 

(8.2) L m : O (1) O (l) ® Zf = M 2 (Z^), 

each C (l)f acts on the left on P^Zp/Ot), and we have £> (l) x /Of ^ ¥ x {Z F /%. 
We conclude that Afjf(tft) = (91) = 0ti-^f (^)i, where 

= {/ : F\Z F /m) W fc (C) :/ | a7 = / for all 7 e O (l)f} ■ 

In this presentation, the Hecke operators act as follows. For a prime p, let 

B(p) id = O (l)f\ {x G hlj 1 : md(xlj^) = p} . 

We then define the linear map T p : .Mf (9T) — > .Mf (9?) on each component by the rule 

(T p )y : («rt), (Sty 

(8-3) /_> £ /| fcT . 

7Ge(p) ii;7 

This is indeed an isomorphism of Hecke modules. For further details, see work of the first 
author [151 Theorem 2] which traces these maps under the assumption that F has narrow 
class number one, but this assumption can be easily removed. 
Put another way, by the decomposition 

h 

B* = \jB*a l d (l)\ 



i=l 



x 



we decompose the set Xq (91) as 

h h 

x ( f(^) = J B x \ J § x /a x = Ufi x \(fi x s i a (i) x )/a x ^□oowrvwr/a 

i=l i=l 

where the last identification is obtained by sending 'jaiU to aiuaj 1 . Thus, analogously 
to (7.17), we get a decomposition 

h h 

(8.4) x B (9t) = [J = □r i \P 1 (z F /*n). 

i=l i=l 

In particular, this gives X(f (1) = Cl(9o(l). From this, we get an identification 
(8-5) ^ 

/ >-K/i)i> 



where we set /(x) = /(««;) after choosing a G Co(l)f such that x = a? • ooj. Again, 
the decomposition (8.5) is analogous to (7.19), and one shows that it is a Hecke module 
isomorphism by arguing similarly. 

Now / G (m) is by definition a map / : B x /O x -»■ Wfc(C) such that /| fe7 = / for 
all 7 G 5 X . Associated to such a map, via the identification (8.4), such a map is uniquely 
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defined by a tuple of maps with fa : 0q(1)*/0* — > Wfc(C) such that / ^7 = / for all 
7 G Tj = C (1) 4 X . In other words, 



M fc B (OT) ^> 0^ o (r j ,Hom(6o(l)rM x ,^(C))) = 0tf°(r,,Coind^ x (1)r W k (C)). 



,Oo(l)f 

i=l i=l 



But then as in (8.2), we have 



H°(r 4 ,CoindJ x (1) > x w k (C)) - (m)*. 



Example 8.6. The real quadratic field F = Q(\/106) has strict class number 2 and class 
number 2. We compute that the space S^l) of Hilbert cusp forms of level 1 and parallel 
weight 2 has dimension 50. It decomposes into four Hecke constituents of dimension 1, six 
of dimension 2, two of dimension 4 and one of dimension 26. The table below contains the 
first few Hecke eigenvalues of the one-dimensional constituents. 



5 


5 


7 


7 


3 


-2 


4 


-1 


-2 


3 


-1 


4 


-3 


2 


4 


-1 


2 


-3 


-1 


4 



p (2,w) (3,w + l) (3,w-l) (5,tw + l) (5,u>-l) (3w + 31) (3w - 31) 

iVp _2 3 3 

a P (fi) -1 -2 3 

a P (/ 2 ) -1 3 -2 

a„(/ 3 ) 1 2 -3 

a„Cf 4 ) 1 -3 2 

The forms /1 and fi (resp. ^3 and /i) are interchanged by the action of Gal(F/Q) (on the 
ideals p). The forms f\ and / 3 (resp. fi and fi) are interchanged by the action of Cl + Zp, 
so these forms are twists via the strict class character of Ga\(F + / F) , where F + denotes the 
strict class field of F. 

Elkies has found a curve E which gives rise to the above data: 

E : y 2 — wxy - 2wy = x 3 + (-2 - 2w)x 2 + (-10809936 + 1049944w) 

+ (-19477819120 + 1891853024w). 

The curve E has j-invariant j(E) = 264235 + 25777u> and has everywhere good reduction. 
We conclude that E is modular using Kisin [421 Theorem 2.2.18] (see also Kisin [43| Theorem 
3.5.5]): we need to verify that 3 is split in F, that E has no CM, and that the representation 
P3 : Gal(F/F) — > GL.2(Z3) has surjective reduction p 3 : Gal(F/F) — > GL^Fs) which is 
solvable hence modular. 

We find that E matches the form /1; so its conjugate by Gal(F/Q) corresponds to fi and 
the quadratic twist of E by the fundamental unit 4005 — 389u> (of norm —1) corresponds to 
/ 3 (and its conjugate to / 4 ). 

The input of our algorithm is a totally real number field F of degree n, a totally definite 
quaternion algebra B with discriminant D, an integral ideal 9T C % F which is coprime with 
D, a weight fceZ" such that K>2 and fc, = kj (mod 2), and a prime p { D. The output 
is then a matrix giving the action of T p in a basis of M^i^ft) = ® i=1 -Mf which is 
independent of p. By computing enough T p and simultaneously diagonalising, one obtains 
all Hecke constituents corresponding to Hilbert newforms of level 91 and weight k. 
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The algorithm starts by finding a maximal order Oq(1), then computes a set of repre- 
sentatives ClCo(l) of the right ideal classes of O (l) whose norms generate C1 + (Z F ) and 
are supported outside Q^l. This part of the algorithm uses work of the second author and 
Kirschmer |41j . it is the most time consuming part but can be seen as a precomputation. 
Next, the algorithm finds a fundamental domain for the action of each Tj on P 1 (Z^/9 r t), and 
computes Mjf (91) as the direct sum of the 

where T x is the stabilizer of x in IV From this, one obtains a basis of M^i^St). Finally, the 
algorithm computes the sets 0(p)jj, and then the block matrices which give the action of T p 
in this basis. We refer to [T5] and [TS] for further details on the implementation. 



9. Indefinite method, arbitrary class number 



In this section, we generalize the indefinite method to arbitrary class number. We carry 
over the notation from Section 7, and now take the quaternion algebra B to be ramified at 
all but one real place. 

In this case, from (T^-frlgl the space X(C) = X B (9t)(C) = B$\(H x B x /<5 X ) is the 



disjoint union of Riemann surfaces indexed by CI Z F . Let {a} be a set of representatives 
for Cl + Ijf and let a G Zj? be such that aZj? n Z^ = a for each a. Then 



(9.1) 



X(C) 



□ r a \n= □ x a (c) 



[a]eCl+(Z F ) 



[a]GCl+(Z F ) 



where O a = aOoT 1 n B and r a = C a x + . 

Therefore, a modular form of weight k and level 9T is a tuple (/„) of functions f a : H 
Wk(C), indexed by [a] E C\ + Z F , such that for all a, we have 

(f a \kl)(z) = f a (z) 

for all 7 G T a and all z G %. In particular, if k — (2, . . . , 2) is parallel weight 2, then (/„) 
corresponds to a tuple of holomorphic 1-forms ((2m)f a (z) dz) a , one for each curve X a (C). 

We compute with this space of functions by relating them to cohomology, and for that 
we must modify the coefficient module. Define the right GL 2 (C) n = GL 2 (C) x GL 2 (C) n_1 - 
module 



V k (C) 



P^K)(C) = P Wl (m x )(C) ® W k {C). 



The group B x acts on 14(C) via the composite splitting B x °->- GL 2 (C) ra given by 7 1— > 
The Eichler-Shimura theorem, combined with the isomorphism (7.19), applied to each 



component X a (C) of X(C) in ( |9.1 ), gives the isomorphism of Hecke modules 
(9.2) 



5f(9t)^0H- 1 (T B) %(C))- 



where + denotes the +l- eigenspace for complex conjugation. 

In the description (9.2), the Hecke operators act on H 1 ^^, 14(C)) in the following way; 



we follow their definition in (7.24). Let p be a prime ideal of lip with p { 2)91. We consider 
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the [b]-summand, and given / G H 1 ^^, Vfc(C)) we will define T p f G H 1 ^^ Vfc(C)), where 
[b] = [p _1 a]. Let I a = aODB and I b = fidnB so that nrd(/ b ) = b and nrd(J ) = a, and let 



e(p)„, b = r 6 \ G ^ n V" 1 : nrd^/^p = (nrd(ro))} 
= T b \ {ro G B* fl /b/" 1 : nrd(ro)a = pb} , 



where r b = O b+ acts by multiplication on the left. Let 7 G r a , so that 7/,, = J a . Then the 
map w 1 — ^ 077 on _B X induces a bijection (of the equivalence classes) of 6(p) 0)b . Therefore, 
for every w G 0(p) 0) b, there exists 5 m G r a and ro 7 G @(p) ,b suc h that 077 = £ ro G7 7 . From 



(7.22) and the Eichler-Shimura theorem, we have 



(9-3) (V)( 7 )= £ 

roee(p) aif) 



One can similarly define the Atkin-Lehner involutions. 

Admittedly, this description is com plic ated, but it can be summarized simply: a Hecke 
operator T p permutes the summands ( |9.2[ ) in accordance with translation by [p] in Cl + Zj7, 
and adjusting for this factor one can principalize as before (when the strict class number was 
1). The resulting Hecke matrices are consequently block matrices. 

We illustrate this with an example; we then give a few more details on the algorithm. 



Example 9.4. Let F = Q(w) be the (totally real) cubic field of prime discriminant 257, with 
w 3 — w 2 — Aw + 3 = 0. Then F has Galois group S3 and 7jp = Z[w]. The field F has class 
number 1 but strict class number 2: the unit (w — l)(w — 2) generates the group Zp + /Zp 2 
of totally positive units modulo squares. 

Let B = ( — j be the quaternion algebra with i 2 = —1, j 2 = w — 1, and ji = —ij. 

Then B has discriminant D = (1) and is ramified at two of the three real places and 
unramified at the place with w h-> 2.19869..., corresponding to : B <^-> M 2 (IR). The 
order 



O = Z F © (w 2 + w - 3)iZ F © ((w 2 + w) - 8i + j) /2Z F © ((w 2 + w - 2)i + ij) /2Z F 



is an Eichler order of level = (w) 2 where N(w) = 3. 

A fundamental domain for the action of F = Loo(0+) on T-L is as follows. 
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The ideals (1) and a = (w + l)Zp represent the classes in the strict class group Cl + Zp. 
The ideal J a = 20 + {{w 2 + w + 2)/2 - Ai + (l/2)j) O has nrd(J a ) = a. The left order of J a 
is Ol{Jo) — Co where 



£> a = Z F © (V - 2w - 3)iZ F © ((V + w)/2 - Ai + (l/2)j) Z F 

© (1/10) ((174w 2 - 3A3w - 348)i + (w 2 - 2w - 2)j + (-w 2 + 2w + 2)ij) Z F . 



A fundamental domain for the action of T a = too(Co+) 011 ^ ^ s as follows. 
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The orders O and O a are not isomorphic since the connecting ideal I a (with left order O a 
and right order O) is not principal. This implies that the groups V and r a are not conjugate 
as subgroups of PSL 2 (IR) but nevertheless are isomorphic as abstract groups: they both have 
signature (1; 2, 2, 2, 2), so that 

r = r a = ( 1 , 1 ',5 1 ,...,5,-.5 2 1 = --- = sl = [ 7 ,7']5i---54 = i). 

In particular, both X^(C) and X a (C) have genus 1, so 

dimiJ 1 (X(C),C) = dimiJ 1 (X (1) (C),C) + dimif 1 (X a (C),C) = 4 = 2dimS 2 (9t). 

We choose a basis of characteristic functions on 7,7' as a basis for if 1 (X( 1 )(C), C) and 
similarly for H\X a (C),C). 

We now compute Hecke operators following the above. Let H = if 1 (X(C), C). We 
compute that complex conjugation acts on H by the matrix 



H\W n 



(-1 





-1 o\ 

1 
0-10 
1 1/ 

Note that W 7 ^ in this case preserves each factor. Now consider the Hecke operator T p where 
p = (2w — 1) and N(p) = 7. Then p represents the nontrivial class in Cl + Z^. We compute 
that 

/ -3 2 \ 
0-2-4 
-4-2 
2 -3 J 
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H I T h 



V 



and restricting we get 

- 2 ). 

Therefore there are two eigenspaces for T p with eigenvalues 4, —4. By contrast, the Hecke 
operator T( 2 ) acts by the scalar matrix 3 on H, preserving each component. 
Continuing in this way, we find the following table of eigenvalues: 



Np 


3 7 


8 


9 


19 25 


37 


41 


43 47 


49 


53 61 


61 


61 


a P (f) 


-1 4 


3 


-4 


-4 -8 


4 


-6 


-8 


4 


12 -8 


2 


4 


ap(g) 


-1 -4 


3 


4 


-4 8 


-4 


-6 


8 


-4 


-12 8 


2 


-4 



Clearly, the form g is the quadratic twist of the form / by the nontrivial character of the 
strict class group Gal(F + /F), where F + is the strict class field of F. Note also that these 
forms do not arise from base change from Q, since a p has three different values for the primes 
p of norm 61. 

We are then led to search for elliptic curves of conductor 9? = (w) 2 , and we find two: 

E f :y 2 + (w 2 + l)xy = x 3 - x 2 + (-36w 2 + 51w - 18)x + (-158w 2 + 557w - 317) 

Eg : y 2 + (w 2 + w + l)xy + y = x 3 + (w 2 - w - l)x 2 + (4w 2 + llw - ll)x + (4w 2 + w - 3) 

Each of these curves have nontrivial Z/2Z-torsion over F, so as above they are modular and 
we match Hecke eigenvalues to find that Ef corresponds to / and E g corresponds to g. 

In this situation, although by the theory of canonical models due to Deligne we know that 
the variety X(C) = X^(C) U X a (C) has a model Xp over F, the curves themselves are 
not defined over F — they are interchanged by the action of Gal(F + /F). Nevertheless, the 
Jacobian of X F is an abelian variety of dimension 2 defined over F which is isogenous to 
Ef x Eg — we characterize in this way isogeny classes, not isomorphism classes. 

As in the case of class number 1, the application of Shapiro's lemma allows us always to 
work with the group associated to a maximal order, as follows. Let Oq(1) D O be a maximal 
order containing O, and for each ideal a, let Oo(l) a = q.Oq{1)q.~ 1 D B be the maximal order 
containing O a , and let r(l)„ = too Further, define 

\4(C) a = Coindg 1)o \4(C) 
for each a. Then Shapiro's lemma implies that 

F 1 (r a ,\4(C)) = /J 1 (r(i) a ,\4(C) a ) 

and so 

s?{vt)^®H\r(i) a ,v k (c) a ) + . 

[a] 

Our algorithm takes as input a totally real field F of degree [F : Q] = n, a quaternion 
algebra B over F split at a unique real place, an ideal 91 C coprime to the discriminant 
D of B, a vector k G (2Z >0 ) n , and a prime p \ SOT, and outputs the matrix of the Hecke 
operator T p acting on the space H = @ b H l (r(l) b , 14(C)b) with respect to some fixed 
basis which does not depend on p. From these matrices, we decompose the space H into 
Hecke-irreducible subspaces by linear algebra. We give a short overview of this algorithm. 

First, some precomputation. We precompute a set of representatives [a] for the strict class 
group Cl + Ijf with each a coprime to p£>9T. For each representative ideal a, precompute 
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a right Oo(l)-ideal I a such that nrd(/ a ) = a and let Oo(l) a = Ol(Io) be the left order of 
I a . Next, we compute for each a a finite presentation for r(l)„ consisting of a (minimal) 
set of generators G a and relations R a together with a solution to the word problem for the 
computed presentation [67J. Then using standard linear algebra techniques, we compute a 
basis for the space W if 1 (r(l) a , V k (C) a ). 



The main issue then is to make the description (9.3) amenable to explicit computation. 
First, compute a splitting l p : Oq{1) M 2 (Zp iP ). Then for each ideal a, perform the 
following steps. 

First, compute the ideal b with ideal class [b] = [p _1 o]. Compute the left ideals 

J a = O a i- l (l t\+OaP 



indexed by the elements a = (x : y) G P 1 (F p ). Then compute the left C b -ideals I b I~ l J a and 
compute totally positive generators w a G O a R -B+ corresponding to 0^w a = I b I~ l J a [%!] . 

Now, for each 7 in a set of generators G a for r a , compute the permutation 7* of P X (F P ) 
[3T3] Algorithm 5.8] and then the elements 5 a = vj a jw~* a for a G P 1 (F p ); write each such 
element 5 a as a word in G b and then apply the formula 

Cr p / B )( 7 )= E AW 1 *- 

The algorithm in its full detail is rather complicated to describe; we refer the reader to 
work of the second author 169 1 for the details. 
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